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Report  Summary 


This  rqxMTt  documoits  work  on  several  subtasks  of  a  nuclear  noonitoring  research 
program  initiated  in  August  1993,  but  terminated  in  December  1993.  For  that  reason, 
several  of  the  reports  included  herein  are  essentially  preliminary  progress  rq>Qrts. 

An  inqxur^t  thrust  of  the  research  program  was  to  develtq)  improvements  to  the 
Center  for  Sdsmic  Studies'  opmtional  system,  the  Intelligent  Monitoring  Syst^  (IMS). 
This  report  includes  a  description  of  a  strategy  for  optimizing  the  station-specific 
parameters  used  for  automatic  processing  of  data  from  imlividual  stations,  and  illustrates 
this  strategy  by  a  case  study  of  its  applicatirm  to  two  stations.  Another  study  reports  on 
the  develi^unent  of  an  in^noved  mean  re^onal  travel  time  curve  ftn*  Pn  i^rcqniate  fm* 
stations  connected  m  the  and  on  a  me^uxi  for  infening  local  correctimis  to  this  mean 
curve.  This  study  also  applies  Deep  Seismic  Sounding  data  from  the  region  to 
demrxistrate  that  phases  commonly  identified  as  Pg  and  Lg  are  actually  PmP  and  SmS, 
and  that  implication  of  this  knowledge  to  epicenter  determination  removes  a  bias  in 
distance  previously  noted  in  IMS  results. 

The  report  includes  a  study  of  an  event  detected  by  the  IMS,  possibly  on  the  Lop 
Nur  nuclear  testing  site  based  on  ^S  Data  alone.  The  inmo'buice  of  data  contributed  by 
a  global  netw<»k  was  demonstrated  by  estabUshing  with  high  confidence  that  the  event 
did  not  occur  at  previously  known  explosion  testing  areas,  and  occurred  at  a  depth 
inconsistent  with  nuclear  testing. 

A  previous  report  described  a  Ground  Truth  Database,  consisting  of  waveforms 
and  carefully  measmed  phase  paramet^  from  events  independently  known  to  be 
earthquakes,  or  known  to  be  explosions  or  explosion-related.  This  rqxnt  up^tes  the 
previous  report,  adding  additional  events  that  may  be  found  to  be  useful  for  event 
discriminatitm  research. 

Hnally,  this  report  describes  work  initiated  to  exploit  data  from  large  mining 
blasts  in  the  Bnmer  Soviet  Union  (  the  Tymyauz  Mine  in  die  Northern  Caucasus)  to 
improve  location  and  discrimination  of  seismic  events.  The  work  was  done  in 
cooperation  with  scientists  of  the  Russian  Institute  of  Physics  of  Ae  Earth  and  Institute 
for  Dynamics  of  the  Geospheres,  Russian  Academy  of  Sciences.  Their  atteiiq)ts  to 
measure  mechanical-electrical-magnetic  phenomena  and  to  provide  precise  data  on  the 
sequence  of  rmplc-fired  explosions  are  included  here.  Additional  work  cm  this  latter 
ttmic,  acconmli^^i^  jointly  with  scientists  of  SMU,  will  be  reported  elsewhere  in  a  report 
on  Contract  F29601-91-D-DB20. 
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'Dining  of  Detection  and  Feature  Extraction  for  the 

GSE'TT-S  IDC 

Anne  Suteau-Henson  and  Ray  Willemann 


Summary 

This  report  describes  a  strategy  devekqped  to  nqiidly  optimize  tiie  parameters 
Ctoning**)  used  for  automated  processing  of  data  from  candidate  seismic  stations  for  tiie 
global  netwoik  proposed  for  GSETT~3,  an  experiment  planned  by  the  Group  d  Scientific 
Experts  of  tiie  Conference  <mi  Disarmament  The  implication  of  tiiis  strategy  to  die  tuning 
of  detecti(»i  and  feature  extraction  is  outlined  here.  The  tuning  of  other  automated  pro¬ 
cesses  (including,  for  example,  phase  idoitification  and  magnitude  estimation)  is  beyond 
the  scope  of  this  report  The  following  topics  are  addressed  here: 


•  objectives  of  the  tuning  process  in  relation  to  the  new  GSE  concqits; 

•  selection  of  an  appropriate  set  of  continuous  waveform  data  to  be  requested  frmn 
each  station  for  initial  tuning  purposes; 

•  automated  si'giud  processing  of  tuning  data,  using  standardized  sets  of  processing 
parameters  C^i^pes”); 

•  analyst  review  of  the  results  of  automated  signal  processing; 

•  optimization  of  the  detector,  and  of  tiie  measuremoit  of  useful  characteristics  of 
detected  signals  CTeature  extraction”),  such  as  polarization  attributes; 

•  tuning  update  as  new  data  become  available. 

The  proposed  strategy  is  illustrated  witii  a  case  study',  the  actual  tuning  of  feature 
extraction  at  two  stations,  as  it  was  performed  for  the  first  prototype  of  the  GSETT-3  biter- 
national  Data  Center  (IDC). 

Introduction 

The  Group  of  Scientific  Experts  (GSE)  of  the  Conference  on  Disarmament  is  plan¬ 
ning  a  third  Technical  Test  of  its  concepts  fOT  a  global  nuclear  monitoring  netwtxk 
(“GSETT-3").  This  network  will  include  three  different  classes  of  stations:  “alpha"  (about 
50),  which  transmit  continuous  waveform  data  to  the  IDC,  where  they  are  us^  fcx  detec- 
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tion  and  location  of  events;  “beta**  (75  to  ISO),  providing  waveform  segments  requested 
by  the  IDC  for  refinement  of  locations  of  legiond  events  formed  fi[x>m  the  alpha  network 
data;  and  **$001010”^  supplying  additional  regional  data  to  further  iminove  die  near  real¬ 
time  bulletin  produced  by  the  alpha-beta  network.  A  small  prototype  netwcnk  has  been 
used  since  May  1993  at  the  International  Data  Center  (IDC)  of  the  Center  for  Seismic 
Studies,  for  an  experiment  referred  to  here  as  **V0”,  Version  0  of  tiie  IDC  Ibstbed  (GSE/ 
US/8S.  1993;  Stead  etal.,  1993). 

Candidate  alpha  and  beta  stations  for  the  GSETT-3  network  need  to  be  evaluated. 
For  each  selected  station,  optimized  parameters  fw  automated  processing  must  be  deter¬ 
mined  (we  refer  to  this  process  as  “tuning”  die  station).  The  advantages  of  determining 
station-specific  processing  parameters  fw  detection,  feature  extraction,  phase  identifica¬ 
tion,  location,  and  magnitude  estimation  have  been  demonstrated  in  previous  work  (e.g., 
Suteau-Henson,  1991;  Riviere-Barbierer  n/.,  1992;  Suteau-Henson  etal.,  1992).  Here,  we 
build  on  this  previous  experience  with  knowledge  acquisition  to  develc^  standardized  and 
automated  tuning  method  that  can  be  used  to  improve  the  IDC  system  performance. 

This  report  describes  the  strategy  developed  for  rrqiidly  evaluating,  selecting  and 
tuning  GSETT-3  candidate  stations.  It  is  based  on  experience  acquired  during  VO.  The  fol¬ 
lowing  requirements  for  tuning  emerge  from  the  new  GSE  concepts: 

•  rapid  knowledge  acquisition:  region-specific  knowledge  for  each  of  the  many 

stations  of  the  proposed  global  network  must  be  rapidly  acquired  and  incmpo- 
rated  into  the  IDC  system; 

•  optmized  automated  processing:  large  volumes  of  data  will  be  processed  near 
real-time  to  produce  event  bulletins,  so  the  best  possible  automated  bulletin  must 
be  obtained  to  reduce  the  human  analysis  load. 

The  strategy  described  here  reflects  our  attempts  at  meeting  these  two  require¬ 
ments.  An  iterative  approach  is  proposed.  First,  an  initial  tuning  is  performed,  uang  an 
appropriate  data  set  of  continuous  waveform  data.  Then,  updates  are  performed,  by  “re¬ 
tuning”  stations  already  used  in  the  experiment,  by  adding  new  data  that  have  become 
available.  This  process  is  repeated  until  a  given  performance  goal  is  achieved  or  no  further 
improvement  can  be  obtained. 

Selecting  Initial  'Rining  Data 

The  first  step  in  tuning  a  candidate  GSETT-3  station  consists  of  requesting  a  set  of 
continuous  waveform  data  from  the  station.  This  is  done  preferably  before  the  station 
starts  being  used  in  the  experiment  (i.e.,  before  the  start  of  data  transmission),  for  the  fol¬ 
lowing  reasons: 
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•  S(n3e  initial  tuning  can  be  petfonxied*  so  that  tuned  parameters  can  be  used  as 
soon  as  the  station  ctmies  on-line; 

•  die  analysts  have  the  opportunity  to  familiarize  themselves  widi  die  station 
before  they  have  to  analyze  data  finom  it  in  near  real-tim^ 

•  die  station  can  be  evaluated. 

The  following  requirements  drive  the  selectimi  of  initial  tuning  data  sets: 

•  for  a  given  station,  the  tuning  data  must  produce  detectitms  (fixun  noise,  real 
seismic  signals,  and  coda)  fcxming  a  sample  rqnesentative  of  the  detecti(»s  that 
will  be  produced  at  this  station  during  the  experiment; 

•  recent  data  are  preferred,  because  they  ate  more  likely  to  reflect  the  conditions 
(instrumentation,  noise,  natural  and  industrial  seismic  activity)  diat  will  prevail 
at  the  station  during  die  experiment; 

•  efiScient  use  must  be  made  of  the  analyst’s  time  during  the  review  of  die  results 
from  automated  processing  of  the  tuning  data,  as  this  analyst  review  is  the  most 
labor-intensive  part  of  the  tuning  process. 

The  first  and  third  requirements  prompt  us,  whenever  feasible,  to  select  qiedfic 
time  periods  for  tuning  purposes.  Uring  bulletin  information  (e.g.,  from  NEIS  and/or 
regional  catalogs),  time  periods  can  be  selected  on  the  basis  of  their  including  a  sizeable 
set  of  rqnesentative  events.  Because  of  the  second  requiiement,  the  most  recent  data 
offering  a  good  event  selection  are  prefetied.  Also,  the  instrumentation  at  the  station  (and 
amfiguration  for  arrays)  during  the  time  period  covered  by  die  tuning  data  must  be  identi¬ 
cal  to  that  in  effect  during  the  actual  experiment. 

Based  on  our  experience,  a  representative  data  set  that  can  be  used  for  initial  tun¬ 
ing  of  detection  and  feature  extraction  should  provide: 


•  at  least  100  noise  detections,  spanning  at  least  a  whole  continuous  24  hours,  (mi  a 
week  day; 

•  at  least  30  detections  of  real  seismic  signals  representing  each  phase  type 
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(regional  P.  regional  S,  teleseisinic). 

The  requirements  for  die  second  categoiy  (real  signals)  are  somewhat  different  for 
alphaaadbeta  stations.  During  the  IDC  VO  experiment,  data  foom  a  particular  teta  station 
were  requested  only  vdien  an  event  had  been  detected  and  located  within  S”  (ff  the  station 
by  die  DX^  based  on  the  a^ha  network  data.  Therefore,  relevant  arrivals  were  local  / 
regional  P  and  S  phases.  Alpha  stations,  however,  can  detect  and  locate  events  from  any 
epicentral  distance,  and  thus  teleseismic  signals  are  also  lelevanL  Because  of  the  density 
of  the  alpha  netwtxk,  emphasis  is  put  on  the  far-regional  range. 

To  be  representative,  the  sample  of  evoits  included  in  the  tuning  data  set  ideally 
should  cover  a  variety  of  distances  (ivithin  the  range  relevant  for  the  station  9pe),  azi¬ 
muths  and  signal-to-noise  ratios  (SNR).  In  practice,  data  sets  biased  toward  a  particular 
epicentral  location  should  be  avoided. 

An  initial  tuning  data  set  can  be  obtained  through  a  request  for  continuous  wave¬ 
form  data  covering  a  time  period  of  a  few  days  to  a  few  weeks.  Although  such  a  data  set 
may  not  always  be  representative  enough,  this  shortcoming  can  be  remedied  by  perform¬ 
ing  tuning  updates  during  the  experiment,  as  described  in  a  later  section. 

Processing  and  Analysis  of 'Dining  Data 


Automated  Signal  Processing 

The  waveform  data  selected  for  tuning  are  run  through  an  automated  signal  pro¬ 
cessor  (JSigPro),  including  a  detector  and  a  post-detection  processor  that  measures  detec¬ 
tion  characteristics  (Bache  et  al.,  1990).  Standard  recipes  are  used  for  a  given  type  of 
station  (such  as  three-component,  broad-band,  sampled  at  20  Hz).  As  detector  optimiza¬ 
tion  is  perfotmed  based  on  results  frcxn  this  process,  the  detector  recipes  should  include  a 
number  of  “beams”  filtered  in  overlapping  narrow  frequency  bands,  with  detectitxi  thresh¬ 
olds  set  low  so  as  to  minimize  the  number  of  missed  teal  arrivals.  This  is  described  in 
mme  detail  in  the  section  on  detector  optimization. 

Analyst  Review 

The  results  fiom  autonuued  signal  processing  are  then  reviewed  by  an  experienced 
seismic  analyst.  The  objectives  of  this  analyst  review  are  as  follows: 


•  perform  final  identification  for  a  representative  sample  of  detections  (noise,  sig¬ 
nal,  and  coda); 

•  evaluate  the  station,  in  terms  of  noise  and  signal  characteristics,  any  specific 
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problems  (noise,  tuning),  seismic  activity. 


In  prKtice.  the  first  t^jective  requires  the  following  stq>s: 

•  analyze  at  least  one  day  of  continuous  data,  identifying  each  detection  (induding 
noise),  and  adding  any  missed  real  sipud; 

•  analyze  time  segments  containing  events  of  interest  for  the  remaining  days,  idea* 
tifying  each  detection,  and  adding  any  missed  real  agnaL' 

To  save  time  and  ease  the  task  of  the  analysts,  it  is  useful  to  noake  available  to  them 
any  relevant  event  bulletins  (NEIS,  regional  catalog).  This  also  helps  insure  the  accuracy 
of  phase  identification  from  analyst  review,  which  is  inqMXtant  because  it  is  die  reference 
against  which  the  performance  of  automated  processing  is  evaluated. 

Optimization 


Detector 

The  first  step  in  tuning  per  se  consists  of  optimizing  the  detector.  A  method  to  per- 
fonn  detector  optimization  was  developed  and  is  summarized  here.  For  every  detection, 
SNR  is  measured  on  every  filtered  detection  **beam**  firom  the  appropriate  stai^ard  recipe, 
udng  the  lowest  threshold  that  might  be  useful  on  this  beam  to  detect  real  signals.  The 
t^timum  thresholds  are  those  that  minimize  the  function: 

n,  n* 

^(0  =  +  ^w^xFt(t,Pf)  +  ^  (fy-lO)VlOOO 

1=1  y=i 

where  t  is  a  vector  of  beam  thresholds,  n  represents  ntnse  detections  and  Pi  real  arrivals, 
is  the  number  of  beams,  np  is  the  number  of  real  phase  types,  and  W}  are  weights  tint 
reflect  the  importance  of  detecting  each  phase  type.  The  functitms  Ff  and  Ft  are  ranqi 
functions  that  increase  when  the  threshold  detect  “false”  arrivals  or  fail  to  detect  “true” 
arrivals,  respectively.  The  purpose  of  the  last  term  is  to  stabilize  die  minimization  algo¬ 
rithm  with  respea  to  thresholds  outside  the  SNR  range  of  the  tuning  data. 

At  three-component  stations,  the  optimization  is  perfomied  using  different  options 
for  the  “beams”  formed  from  the  vertical,  north-south  and  east-west  channels  (Z,  N,  E, 
respectively),  and  the  rectified  sum  of  the  horizontal  channels  (H).  For  example,  the  fol¬ 
lowing  beams  can  be  used: 
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•  Z,  N,  E,  H,  with  duesholds  set  independently; 

•  Z.  N,  E,H,  constraining  die  thresholds  for  N  and  E  to  be  equal; 

•  Z,  H,  with  independent  thresholds. 

Qeaiiy,  not  all  of  these  options  have  die  advantage  ci  avrading  azhnuthal  bias. 
Directional  preferences,  however,  may  be  desirable  for  some  stations,  to  take  advantage  of 
any  knowledge  about  specific  directions  for  noise  or  seismic  sources. 

Feature  Extraction 

In  order  to  optimiae  automated  phase  ideatification,  it  is  important  to  first  pexfixm 
the  best  possible  measurement  of  characteristics  Creature  extractioa*^  for  eadi  detection. 
A  feature  extractimi  process  is  considered  as  (^timized  when  it  results  in  the  best  dassifir 
cation  of  detections  of  different  types  (ie..  noise,  regional  P,  regi<»al  S,  and,  ice  alpha  sta¬ 
tions,  teleseismic),  the  identification  given  by  the  analyst  being  used  as  reference. 

Fdr  an  array  station,  the  most  useful  feature  extraction  process  for  phase  identifica¬ 
tion  is /-^  analysis  (ot  any  other  process  deriving  the  slowness  vector  from  the  2-D  wave- 
field).  For  a  three-conqionent  station,  the  equivalteit  |»ocess  is  polarizatitm  analysis. 

The  optimization  ccuisists  of  explming  the  parameter  sqiace  for  die  recipe  used,  and 
sdecting  the  set  of  parameters  that  results  in  the  best  classtficatimi  performance.  In  prac¬ 
tice,  our  experience  with  data  from  a  variety  of  stations  makes  it  possible  to  place  reason¬ 
able  limits  on  the  range  of  die  various  parameters  explored.  To  tqitimize  polarization 
analysis,  for  example,  a  set  of  45  standard  redpes  has  bwn  developed  for  data  sanqiled  at 
up  to  80  Hz  (see  case  study  bdow). 

The  data  set  of  detections  used  for  optimization  of  feature  extraction  is  a  subset  of 
that  produced  by  the  standard  detector  us^  fen’  detector  optimization,  ft  includes  only 
those  noise  and  coda  detections  that  passed  the  detection  test  of  the  optimized  detector;  as 
they  are  those  relevant  for  phase  identificatiem. 

In  practice,  this  <q>timizati<m  requires  the  following  st^s: 

•  re-run  the  post-detection  processor  (such  as  polarization  or  f~k  analysis)  off-line 
for  each  detection  in  die  selected  set  and  for  each  of  die  standard  recipes; 

•  for  each  recipe,  automatically  classify  each  detection  using,  for  example,  classi¬ 
cal  discriminant  analysis  (Suteau-HenSon,  1991); 

•  select  the  recipe  with  the  highest  classification  performance. 

Following  tuning  of  feature  extraction,  tuning  of  automated  phase  identificatimi 
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can  be  peifomied  using  a  set  of  detections  (and  their  characterisdcs)  derived  after  optimi¬ 
zation  of  detection  and  feature  extraction.  For  example,  tuning  of  initial  phase  identifica¬ 
tion  at  a  station  of  the  three-component  or  small-array  type  consists  of  training  a  neural 
network  (Patnaik  et  al.,  1992;  Seieno  and  Pamaik,  1993).  As  this  q)timization  process 
requires  a  larger  set  of  detections  in  each  phase  type  than  the  previous  two,  it  is  performed 
when  enough  data  (from  the  initial  tuning  set  and  from  the  experiment  itself)  have  bectxne 
available  for  a  given  station.  The  process  of  re-tuning  after  adding  new  data  to  die  initial 
data  set  to  increase  the  sample  size  is  described  in  more  detail  in  the  next  section.  Further 
description  of  the  tuning  of  phase  identification  is  beyond  the  scope  of  this  work. 

'Hining  Update 

In  the  previous  sections,  we  addressed  the  tuning  of  a  station  using  an  initial  data 
set  As  more  data  from  this  station  become  available  during  the  experiment  itself  (ccmtinu- 
ous  data  from  alpha  stations  and  event  segments  from  beta  stations),  it  is  advantageous  to 
update  the  tuning  by  adding  such  data  to  the  initial  set,  because: 


•  this  increases  the  size  of  the  tuning  data  set,  also  making  it  possible  to  tune  pro¬ 
cesses  requiring  larger  data  sets,  such  as  phase  identification; 

•  this  increases  the  quality  of  the  tuning  data  set,  as  data  from  die  actual  experi¬ 
ment  are  more  relevant; 

•  the  new  data  are  already  analyzed  as  part  of  routine  IDC  operations. 

In  particular,  re-tuning  is  necessary  when  upgrades  have  been  made  to  a  station.  In 
order  to  avoid  inconsistencies,  automated  processes  should  be  re-tuned  in  the  same  order 
as  that  in  which  these  processes  occur  in  the  operational  data  processing  pipeline.  For 
example,  the  tuning  of  initial  phase  identification  (e.g.,  training  of  station-spe^c  neural 
networks)  should  be  performed  using  data  derived  from  running  an  optimized  detectin’  and 
an  optimized  feature-extraction  process. 

A  Case  Study:  Polarization  Hining  of  GRSN  Stations  MOX  and  BRG 

The  proposed  tuning  method  is  illustrated  with  the  optimization  of  polarization 
analysis  at  two  three-component  stations  from  the  German  Regional  Seismic  Network, 
“GRSN”  (GSE/Federal  Republic  of  'Gcrmany/39,  1991).  These  stations,  MOX  and  BRG, 
are  located  about  2°  north  of  the  GERESS  array,  in  an  area  with  a  high  level  of  industrial 
seismic  activity.  During  VO,  they  have  been  used  as  beta  stations,  the  nearest  alpha  station 
being  GERESS.  For  the  purpose  of  initial  tuning,  two  and  three  days  of  continuous  data 
were  obtained  from  MOX  and  BRG,  respectively.  The  tuning  was  later  updated  by  adding 
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events  from  fifteen  days  of  data  recorded  during  VO,  within  the  period  from  May  18  to 
June  20, 1993.  Ingure  1  shows  the  locations  of  the  events  recmded  at  MOX  and  BRG  dur¬ 
ing  that  period  of  the  experiment 

Polarization  analysis  was  optimized  by  comparing  the  performance  foe  classifying 
noise  detections,  region^  P  and  regional  S,  of  45  recipes,  sparuiing  a  wide  range  of  pro¬ 
cessing  parameters.  The  most  sensitive  parameter  was  found  to  be  the  finequency  barxl 
used  in  the  analysis.  Bands  including  frequencies  up  to  2, 4, 8, 16,  and  32  Hz  were  tried. 
The  high  sample  rate  (80  Hz)  made  it  possible  to  investigate  the  usefulness  of  high  fre¬ 
quencies  (16-32  Hz).  Note  that  only  results  firom  polarization  analysis  were  used  in  the 
classification,  and  that  including  other  information  (such  as  context)  would  irxqnove  the 
resulting  performance. 

The  classification  performance  of  the  45  recipes  was  compared  for  MOX  before 
and  after  the  tuning  update  (Figure  2 ).  Hgure  3  compares  the  tuning  results  fm^  MOX  and 
BRG  after  the  update.  Table  1  gives  the  number  of  detections  used  for  optimization  in 

Table  1:  Number  of  Detections  for  Polarization  Timing  at  MOX  and  BRG 


Station 

Noise 

Regional  P 

Regiorud  S 

MOX 

(initial  tuning) 

41» 

21 

44 

MOX' 

(initial  tuning  +  VO) 

81 

86 

80 

BRG 

(initial  tuning  VO) 

80 

59 

88 

a.  half  of  the  available  noise  detectkxis  were  used 


each  case.  At  MOX,  the  update  confirms  the  tuning  results  based  on  the  initial  set  (duu 
was  somewhat  deficient  in  regional  P  phases),  with  a  better  performance  and  a  higher  con¬ 
fidence  (Figure  2).  Figure  3  shows  that  at  BRG,  classification  performance  is  optimized 
when  only  frequencies  up  to  4  Hz  are  included  in  the  analysis.  Such  importance  of  low 
frequencies  at  this  station  (the  teleseismic  band)  is  in  agreement  with  the  observation  drat 
it  is  the  GRSN  station  contributing  the  most  to  the  ISC  (Manfred  Hunger,  personal  com¬ 
munication).  At  MOX,  the  use  of  frequencies  above  8  Hz  does  not  improve  the  classifica¬ 
tion  performance. 

A  second  update  was  performed  fbr  BRG  by  adding  data  obtained  during  the 
experiment  within  the  period  from  June  21  to  October  27, 1993.  This  update  confirmed  the 
results  of  the  previous  one,  which  indicates  that,  in  this  case,  no  further  improvement  was 
obtained  by  re-tuning  again. 
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Figure  1:  Locations  of  regional  events  recorded  at  beta  stations  MOX  (squares)  and 
BRG  (crosses)  during  15  data  days  of  VO.  Triangles  denote  MOX  and  BRG, 
and  the  circle  indicates  the  GERESS  array,  used  as  an  alpha  station. 
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CORRECT  PHASE  ID  (%)  CORRECT  PHASE  ID  (%) 


Figure  2:  Boxplots  describing  the  distributions  of  the  percentage  of  correct  phase  iden¬ 

tification  at  station  MOX  for  polarization  recipes  grouped  by  high  cutoff  of 
the  frequency  band  used  in  the  analysis.  The  top  plot  is  for  the  initial  tuning 
set.  The  bottom  plot  is  for  a  data  set  including  the  initial  tuning  set  and  IS 
days  of  data  obtained  during  VO. 
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CORRECT  PHASE  ID  (%)  CORRECT  PHASE  ID  (%) 


Figure  3:  Boxplots  describing  the  distributions  of  the  p«centage  of  correct  phase  iden¬ 
tification  for  polarization  recipes  grouped  by  high  cutoff  of  tiie  frequenc> 
band  used  in  the  analysis.  Both  initial  tuning  data  and  VO  data  were  used,  foi 
station  MOX  (top),  and  BRG  (bottom). 
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Evaluation  and  polarization  tuning  were  also  perfonned  for  other  three-component 
statitms  used  as  beta  stations  during  VO,  as  desaibed  in  Suteau-Henstm  et  al.,  1^3.  In 
particular,  ten  GRSN  stations  were  evaluated,  and  six  were  selected,  including  MOX  and 
BRG.  At  five  GRSN  stations,  a  large  enough  data  set  was  available  to  tune  polarization.  A 
cmnparison  of  classification  performance  as  a  function  of  the  frequency  band  used  in  die 
analysis  shows  that,  for  aU  stations  except  BRG,  the  best  results  are  obtained  when  using 
frequencies  up  tt>  8Hz. 

Conclusions  and  Reconunendations 

The  prc^msed  strategy  for  tuning  candidate  stations  for  GSETT-3  has  the  advantage 
of  helping  meet  the  requirements  for  rtqiid  knowledge  acquisition  and  optimized  auto- 
mamd  processing  of  large  volumes  of  di^  An  iterative  and  mostly  autoinated  approach 
makes  it  possible  to  start  tuning  a  new  statim  befcxe  it  is  incmponued  into  die  network, 
and  progressively  inqirove  this  tuning  as  more  data  becmne  available  during  the  experi¬ 
ment  Tuning  updhues  can  also  be  performed  by  using  new  recipes  or  new  algorithms,  with 
the  aim  of  further  iziqnoving  the  autmnated  processing.  In  practice,  we  expect  that  at  srane 
point  the  performance  inqirovement  obtained  by  re-tuning  a  given  station  will  become 
negligible.  Then,  no  re-tuning  will  be  needed  unless  some  changing  conditions  at  the  sta¬ 
tion  call  for  further  updates. 

The  prt^iosed  rqipioach  was  described  in  m(»e  detail  for  the  case  of  optimization 
of  detection  and  feature  extraction.  It  is,  however,  a  general  approach  that  can  be  applied 
to  the  tuning  of  other  autcmiated  processes  as  well.  Also,  “metrics”  can  be  develop^  to 
objectively  quantify  the  inqnovement  achieved  after  tuning.  This  is  the  object  of  current 
and  future  wmk  aiming  at  improving  the  IDC  perfcumance. 
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A  Note  on  an  Event  in  the 
Xinjiang  Province,  China  on  April  14, 1993 

Hans  Israelsson 


Introduction 

On  April  14  1993  IMS  detected  and  located  a  small  event  (mb=4.3)  with  an  enor 
ellipse  that  covered  the  Lop  Nur  Nuclear  We^KMis  Test  Base  in  the  ^jiang  Brovince  oSl 
China.  Because  of  the  uncertainty  of  the  location  and  the  proximity  of  the  Lop  Nur  test 
base  to  active  seismic  regions  the  nature  of  the  source  caimot  readily  be  established  on  die 
location  alone.  Furthermore,  the  waveforms  at  the  IMS  arrays  were  more  complex  than 
those  of  previously  observed  nuclear  explosions  at  Lop  Nur  (Bratt,  1993). 

In  this  note  we  have  attempted  to  reduce  the  uncertainty  of  the  IMS  location  and 
ounpile  data  to  facilitate  identificadon  of  the  source.  In  doing  so  we  were  drawing  uptm 
data  recorded  at  additional  stations  that  were  requested  from  Natitmal  Data  Centers  in  sev¬ 
eral  countries  and  retrieved  from  the  IRIS/DMC. 

Data 


The  map  in  Figure  4  shows  the  95%  confidence  ellipse  of  the  epicenter  estimated 
from  arrival  times,  azimuths,  and  slowness  values  of  primary  and  secondary  phases  at  four 
of  the  IMS  arrays,  APAO,  ARAO,  FIAO,  and  NRAO  (Table  2)  as  available  in  the  IMS  data 
base  at  the  Center  for  Seismic  Studies.  The  area  of  the  ellipse  is  about  900,  OCX)  km^. 

The  map  also  shows  the  four  IMS  arrays  and  other  stations  which  recorded  tiie 
event  These  include  the  arrays  ASAR  and  WRA  in  Australia  and  YKA  in  (Canada,  all  at 
teleseismic  distances,  and  the  station  WMQ  in  die  Xinjiang  Province  of  China,  at  a  local 
distance  fiom  the  preliminary  IMS  epicentor  estimate. 

Waveforms 

Hgure  5  shows  the  short  period  beams  recorded  at  the  IMS  arrays.  The  beam  for 
NRAO  is  unfiltered  whereas  the  beams  of  the  other  arrays  were  filtered  between  l.S-4.0 
Hz.  The  arrival  times  of  the  first  and  secondary  phases  as  reported  in  the  IMS  data  base, 
are  marked  in  the  record  sections.  The  secondary  phases,  detected  on  the  beams  at  all  four 
arrays,  albeit  with  varying  SNR,  are  delayed  about  11  s  relative  to  the  first  arrival  time. 
The  f-k  spectra  of  the  primary  and  secondary  phases  were  in  general  agreement  suggesting 
that  the  two  phases  originated  from  the  same  source  region. 

SNR  and  dominant  frequency  of  the  primary  and  secondary  phases  vary  signifi¬ 
cantly  among  the  arrays.  The  beam  at  NRAO  had  the  largest  SNR  and  ^e  initial  P  signal  at 
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BEAMS  AT  IMS  ARRAYS  FOR  93/0<4^/1>4  event 


1 0  20  30  40  SO 

TIIME  <S) 


HgureS:  Short  period  phase  beams  of  the  event  on  ^ail  14, 1993,  at  the  four  IMS 

airays,  APAO,  ARAO,  HAO,  and  NRAO.  The  beam  at  NRAO  is  unfiltered 
wheteas  the  beams  at  odm  arrays  were  band  pass  filtered  between  1^ 
4.0  Hz. 
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'Hible  2:  Reported  arrival  data 


Station 

Phase 

Hme 

SD 

Azimuth 

SD 

SlOWTL 

SD 

APAO 

P 

8:38:23.7 

2.3 

93.2 

16.8 

APAO 

Tx 

8:38:35.3 

3.8 

100.7 

6.4 

... 

— 

FIAO 

P 

8:38:43.4 

84.1 

22.6 

6.86 

5.72 

ARAO 

P 

8:38:49.3 

1.0 

87.5 

14.2 

8.70 

5.06 

NRAO 

P 

8:39:39.4 

1.0 

77.1 

13.8 

6.22 

3.53 

NRAO 

Tx 

8:39:52.1 

1.0 

76.7 

13.6 

— 

— 

YKA 

P 

8:42:38.3 

1.0 

WMQ 

P 

8:31:30.0 

2.0 

— 

WMQ 

S 

8:31:44.7 

2.0 

- 

WRA 

P 

8:42:52J 

1.0 

— 

WRA 

pP 

8:43:01.2 

1.0 

— 

ASAR 

P 

1.0 

— 

ASAR 

pP 

8:43:16J 

1.0 

— 

— 

NRAO  had  a  strikingly  lower  dominant  frequency  than  those  of  the  other  arrays. 


Waveform  segments  of  all  array  elements  and  beams  and  f-k  q)ectra  were  received 
via  FAX  for  the  Australian  arrays  ASAR  and  WRA  (Muirliead  and  Penders,  private  com¬ 
munication,  1993),  which  are  reproduced  in  Hgure  6.  Again  clear  secondary  phases  are 
developed,  however,  now  delay^  with  <xily  abmit  8  s  after  the  initial  P.  As  for  die  IMS 
arrays  die  f-k  spectra  of  the  primary  and  secmidary  phases  were  in  general  agreement 

Three  conqionent  recordings  at  the  Chinese  station  WMQ  retrieved  from  die  IRIS 
DMC  are  ^own  in  Figure  7  .  The  large  amplitude  S  wave  dominates  the  recordings.  The 
initial  P  wave  had  a  rectilinearity  close  to  1.0  widi  back  azimuth  (about  190^  and  inclina¬ 
tion  angle  (20*^  well  defined.  The  S  phase  is,  cm  the  vertical  and  north  component  pre¬ 
ceded  with  about  1  s  by  a  phase  with  high  rectilinearity  and  also  well  defined  back 
azimuth  and  angle  of  incUnation  similar  to  that  of  the  initial  phase. 

Other  Data 

A  first  arrival  time  at  the  array  YKA  in  (Canada  was  provided  along  with  locaticm 
and  magnitude  (North,  private  communication,  1993).  The  IMS  array  GEAO  did  not  detect 
the  event  automatically  and  signals  that  could  be  positively  associated  with  the  event 
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Hgure6: 
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.Single  element  traces  and  beams  filtesred  and  unfilteied  at  die  Austra¬ 
lian  array  ASAR  of  the  event  on  April  14, 1993.  The  f-k  plectrum  for 
the  time  window  indicated  on  the  filteied  beam  is  shown  to  die  rig^t 
in  the  figure. 


19 


REOORDINGS  AX  WMO  OF  03/04/1  >4-  EVEMX 


Hgnre  7:  Three  component  teoad  band  reoofdings  at  the  Chinese  station  WMQ  of 

the  event  on  April  14, 1993. 


20 


could  neither  be  found  in  a  subsequent  analysis  (Henger,  private  communicadon.  1993). 
At  the  IRIS  DMC  there  were  no  data  for  a  dme  window  of  expected  arrival  times  for  P  or 
S  phases  at  any  Chinese  stations  other  than  WMQ.  A  request  was  also  made  for  data  at  the 
station  NIL  in  Pakistan,  which  recorded  the  Lop  Nur  explosion  of  92/09/2S  with  a  veiy 
large  signal-to-noise,  but  so  far  no  response  has  been  received. 

Location 

All  locations  were  made  with  the  LocSAT  program  that  is  employed  by  the  IMS. 
In  addition  to  the  epicenter  the  program  estimates  a  confidence  ellipse  tluu  dq)ends  partly 
on  the  a  priori  or  assumed  standard  deviations  of  the  errors  of  measured  arrival  times, 
back  azimuths  and  slowness  values.  The  standard  deviations  in  the  CSS  data  base  were 
used  for  the  measurements  at  the  IMS  arrays.  The  arrival  times  at  the  other  stations  were 
assigned  values  listed  in  Table  2.  As  the  a  priori  standard  errors  of  back  azimuths  and 
slownessAncidence  angle  for  ASAR  and  and  for  WMQ  could  not  be  readily  esti¬ 
mated,  only  the  arrival  times  at  these  stations  were  used  in  the  event  locations. 

The  95%  confidence  ellipse  in  Rgure  4  of  the  initial  IMS  solution  is  based  only  on 
data  at  four  IMS  arra ys  including  arrival  times,  slownesses  and  back  azimuths  for  primary 
and  secondary  phases.  No  convergence  was  obtained  if  the  IMS  data  were  limited  to  first 
arrival  times.  Figure  S  shows  how  the  confidence  ellipse  is  reduced  -  the  length  of  the 
major  axis  goes  from  1200  to  35  km  -  as  data  at  the  other  stafions  are  added  in  die  order 
(1)  YKA,  (2)  WMQ,  and  (3)  ASAR/WRA.  The  standard  deviation  of  the  arrival  time 
residuals  for  the  solution  based  on  all  available  station  data  was  0.62  s. 

Phases  P  and  S-  instead  of  Pg  and  Sg  -  were  used  for  WMQ  in  the  location  as  die 
IMS  travel  time  tables  for  Sg  did  not  extend  to  depth  below  5  km.  Solutions  with  Pg  and  S 
yielded  very  similar  results.  Figure  6  compares  the  velocity  model  used  by  LocSAT  fix’ 
local  and  regional  stations  with  that  obtain^  by  Barker  (1992)  for  SSE  fiom  WMQ. 

The  area  of  the  95%  error  ellipse  for  the  location  based  on  all  station  data  is  about 
1300  km^,  and  the  standard  errors  in  latitude  and  longitude  were  7  and  13  km  respectively. 
These  standard  errors  are  comparable  to  those  reported  by  NEIS  for  its  epicenter  determi¬ 
nations  of  earthquakes  in  the  area  with  magnitude,  mb,  less  than  5. 

Data  With  Bearing  on  the  Nature  of  the  Source 

In  order  to  facilitate  ascertaining  the  nature  of  the  source  of  the  93/04/14  event,  we 
have  compiled  data  according  to  a  simple  “screening”  principle  comnoonly  used  for  event 
identification:  (1)  epicenter,  (2)  focal  depth,  and  (3)  characteristics  of  recorded  wave¬ 
forms. 
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5  6  7  8 

VELOCITY  (km/s) 


Figure  9:  Comparisons  of  veloci^  models  used  by  the  location  program  LocSAT 
and  estimated  by  Baiker  (1992)  foe  a  profile  SSE  the  stati(»  WMQ. 
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Epicenter 

We  conqiare  the  estimated  epicenter  and  associated  uncertainty  with  die  locations 
of  the  Lc^  Nur  nuclear  weapons  test  base  and  ndghboring  seismic  ztxies. 

JED  of  Lop  Nur  explosions 

I^om  the  mi^  in  Hgure  8  it  is  clear  thid  the  ctmfidence  ellipse  of  locatioo  for  the 
93/04/14  is  well  out^kle  the  area  planned  by  die  epicenters  of  previous  nuclear  eiqilosions 
at  the  Lqp  Nur  base,  which  covers  an  area  of  more  than  100, 000  km?.  The  expl(»k»  epi¬ 
centers  in  Hgure  8  were  obtained  with  the  JED  method  (Marshall  and  Jones,  private  com- 
municadon,  1993).The  in  Hgure  10  zomns  in  on  the  explosion  qiicenters  showing 
that  the  atmospheric  and  underground  tests  are  separated  with  almost  100  km.  The  under¬ 
ground  tests  group  in  three  areas  separated  by  about  40  kriL  The  explosions  of  the  most 
eastern  area  (6  events)  are  compared  in  Hgure  10  with  the  locadons  of  two  explosions 
(denoted  Yongtai  in  the  Figure  10).  both  supposed  to  have  occurred  prior  to  1987  as 
described  by  Matzko  (1992)  citing  a  Chinese  publication.  The  JED  of  the  events  in  Hgure 
10  were  constrained  to  the  qiicenter  of  evmt  87/06/05  marked  in  Hgure  10.  Frtmi  this  it 
appears  that  the  events  dated  84/1Q/03  and  83/10/06  in  Hgure  10  could  well  conrespcMid  to 
the  ones  specified  by  the  Chinese,  as  they  ate  ofiEset  only  a  few  Iotl  Therefore,  even  if  die 
uncertainties  of  the  JED  solutions  are  not  available,  the  solutions  appear  to  accurate  at 
least  to  within  a  few  km. 

Seismicity  and  Faults  in  Eastern  Tien  Shan 

The  Lop  Nur  test  base  in  near  the  Hen  Shan  mountain  range,  which  is  a  seismi- 
cally  active.  Important  thrust  faulting  occurs  in  the  north  eastern  Hen  Shan  (Tiqiponier  and 
Molnar,  1979).  Mechanisms  for  large  earthquakes  near  Uninochi  -  where  the  station  WMQ 
is  located  -  show  that  faulting  consistent  with  a  P  axis  roughly  perpendicular  to  the  recent 
E-W  fold  trends.  The  maps  in  Figure  11  and  12  show  the  geographical  distributicm  of  the 
seismicity  since  1970  as  reported  by  USOS  for  an  area  -f-  3^  and  -k  6°  in  latitude  and  Icm- 
gitude  fixim  die  center  of  the  Lop  Wur  test  base.  Major  active  faults  are  also  overlain  in 
Hgure  11  (after  T^ponier  and  Molnar,  1979).  Hie  magnitude  distribution  of  the  events  in 
Hgure  13  shows  foat,  in  average,  diete  were  since  1970  a  few  events/year  larger  than 
mb»4.S.  Large  histtnical  earthquakes  have  also  been  rqxirted  for  this  part  of  the  Hen 
Shan  mountains. 

For  comparison  the  confidence  ellipse  of  the  93/04/14  event  is  also  marked 
together  with  the  earthquakes  located  by  NOS  in  Figure  11  and  is  covers  a  major  fault 
Most  of  the  NEIS  locations  do,  however,  not  line  up  along  the  major  faults  marked  in  Hg¬ 
ure  11.  As  the  standard  errors  of  the  NEIS  epicenter  determinations  are  conqiarable  to 
those  obtained  for  the  93/04/14  here  the  symbol  sizes  in  Hgure  11,  which  ate  proportional 
to  event  magnitude,  should  be  much  larger  to  reflect  location  uncertainty.  It  therefixe 
appears  that  the  epicenter  estimate  of  the  93/04/14  is  entirely  consistent  with  the  prevail¬ 
ing  seismic  activity  in  the  Eastern  Hen  Shan. 
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UCD  or  i.OR  t^UR  eXRI-OSIOtMS 


OOlsJS^RAIlMlfMO  UCO  UOI»  MOR  e:XRI.ORIOl^9 


Figure  10:  (a)  distribution  of  explosion  epicenters  at  the  Lop  Nor  test  base,  under¬ 

ground  and  atmospheric  tests  are  indicated  with  crosses  and  triangles 
respectively;  (b)  relative  location  of  events  used  to  constrain  the  JED  of  the 
explosions. 
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EARTHQUAKES  AND  FAULTS 


PROJECTION:  LJUJSCRT.  CENTER:  «2.a8  07.00  0.0 

WINDOW  ANCt.ES  2.7010  2.7010  2.7010  2.7010 


Hgnre  11:  The  qinoemral  distribution  of  the  seismicity  since  1970  as  lepoited  by 
uses  in  the  Lop  Nor  and  nei^boring  areas.  Major  active  faults  are 
defined  by  Tey)ponier  and  Molnar  (1979)  are  also  indicated  togetho-  with 
the  95%  confidence  ellipse  of  tiie  event  on  April  14, 1993. 
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DISTRIBUTION  OF  mb  FOR  EARTHQUAKES  AROUND  LOP  NUF 


Hguxe  13:  The  magnitude,  inb(USGS),  distribution  of  earthquakes  iqKXted  by 
uses  since  1970  in  the  Ixtp  Nor  and  neighboring  axeas. 
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The  map  in  Hguie  12  provided  by  Fielding  (private  communicaticw,  1993)  shows 
that  the  location  of  the  93/04/14  is  in  quite  a  mountainous  part  of  the  Hen  Shan. 

Interestingly,  some  of  the  previous  earthquake  epicenters  fall  inside  oc  very  close 
to  the  Lop  Nur  test  base  (see  Table  3). 


Table  3:  Earthquakes  near  Lop  Nur  Test  Base 


Date 

Hme 

Lat(N) 

Lon(E) 

Depth 

mb 

76/03/20 

04:34:00 

41.7640 

88.6700 

.  8 

5.00 

87/12/22 

00:16:39 

41.3600 

89.6380 

21 

5.90 

88/11/15 

16:56:46 

42.0180 

89.2950 

33 

5.00 

90/01/21 

07:53:32 

41.5340 

88.7280 

33 

4.60 

Focal  Depth 

It  is  well  known  that  an  accurate  detmnination  of  focal  depth  by  itself  can  posi¬ 
tively  identify  a  seismic  event  The  secondary  phases  at  the  ASAR  and  WRA  and  the  IMS 
arrays  fit  expected  arrival  times  of  pP  and  sP  respectively  which  constrain  the  focal  dq>th 
to  between  2S-30  km,  which  is  compatible  with  the  P  and  S  arrival  times  at  WMQ  and  first 
arrival  times  at  all.stations.  Furthermore  the  similarity  of  the  f-k  spectra  for  primary  and 
secondary  phases  at  the  arrays  lends  support  to  this  interpretation  of  the  secondary  phases 
at  the  arrays  as  depth  phases.  The  estimated  depth,  27-4-2  km,  is  close  to  the  average 
depth  of  earthquakes  in  eastern  Hen  Shan  (see  Figure  14).  Omitting  depth  phases  the  esti¬ 
mated  depth  became  21  +-20  km  and  the  standard  deviations  of  the  residues  is  similar  to 
that  of  the  residuals  for  the  depth  constrained  by  the  depth  phases. 

Characteristics  of  recorded  waveforms 

In  order  to  compare  the  waveforms  recorded  for  the  93/04/14  with  previous  events 
the  data  base  at  CSS  and  the  IRIS/DMC  were  searched  for  events  at  at  close  to  Lop  Nur 
recorded  by  the  IMS  arrays  and  WMQ.  Hgure  13  shows  the  qiicenter  of  the  events  which 
are  summarized  in  Table  4  relative  to  the  location  of  the  93/04/14  event 

Data  at  NRAQ 

Figure  15  compares  waveforms  at  NRAO  (single  sensor)  for  the  93/04/14  event 
with  those  of  a  presumed  underground  nuclear  explosion  (92/09/25)  and  three  other  earth¬ 
quakes.  The  explosion  record  is  on  the  whole  simpler  and  has  a  higher  fiequency  content 
that  do  the  earthquake  records  and  that  of  the  93/04/14  event,  although  these  records  also 
have  clear  high  frequency  components  over  riding  the  dominant  low  frequency  ampli¬ 
tudes.  The  apparently  reduced  low  frequency  components  of  the  explosion  record  agrees 
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DISTRIBUTION  OF  DEPTH  FOR  EARTHQUAKES  AROUND  LOP  NUR 


Figure  14;  The  distribution  of  focal  depths  reported  by  USGS  for  earthquakes  in 
eastern  Tien  Shan. 
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EARTHQUAKES  AND  92/09/25  EXPLOSION  RECORDED  AT  NRAO 


PROJCCnON:  lAMBERT.  CENTER:  42.01  87.S7  0.0 

WINDOW  ANCLES  1.8010  1.8010  1.8010  1.8010 

Hgure  IS:  Epicenters  of  some  events  at  and  near  the  Lop  Nor  test  base  recorded  at 
die  IMS  array  NRAO. 
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Tible  4:  Earthquakes  with  data  for  waveform  oompariaoiis 


Date 

Time 

Lat(N) 

Lon(E) 

Dq>fo 

mb 

91/06/06 

08d)2.*08 

42.705 

87.221 

32 

5.10 

91/17/91 

18:53:22 

43.141 

87.968 

22 

4.80 

92/06/08 

09:20:55 

43.598 

88.277 

30 

4.20 

92/11/27 

im 

41.978 

89.283 

14 

5.30 

with  the  results  of  Ihylor  and  Marshall  (1991)  and  *rsvang  et  aL  (1993).  The  lack  of  low 
frequency  components  of  the  explosion  rectxd  and  its  relative  sinq>licity  are  in  quantita¬ 
tive  agreement  with  the  results  of  Arara  and  Basu  (1984)  who  c(M:iq>ared  ^ectral  and  com¬ 
plexity  parameters  from  recordings  at  the  array  GBA  in  India  from  earthquakes  in 
Xinjiang  and  Lop  Nur  explosions. 

The  difference  in  frequency  content  is  largest  between  the  explosion  record  and  the 
one  for  the  earthquake  closest  to  the  explosion  (within  about  50  km).  The  wavefonn  of  the 
93A)4/14  event  resembles  die  most  that  of  die  closest  earthquake  (91/06/06). 

DataatWMO 

Recordings  at  WMQ  could  only  be  found  for  earthquakes  in  IRIS/DMC  The 
waveforms  for  two  earthquakes  are  compared  with  those  of  the  93/04/14  event  in  Hgure 
16.  The  two  earthquakes  are  at  somewhat  closer  distances,  but  the  associated  waveforms 
appear  quite  similar  but  widi  regard  to  frequency  content  and  reladve  P  and  S  anqilitudes. 
Indeed  foe  S/P  amplitude  ratio  for  foe  93/04/14  on  foe  vertical  componoit  is  stxnewhat 
larger  than  for  foe  two  earthquakes  and  has  a  high  value  typical  of  earthquakes  in  other 
regions. 

Concluding  Remarks 

The  analyzed  here  for  foe  93/04/14  event  all  point  to  an  earthquake  source,  with 
regard  to  location,  focal  depth  as  well  as  recorded  waveform  characteristics. 

The  collected  data  also  illustrate  foe  value  of  high  sensitive  arrays  for  detecticm  of 
small  events.  The  availability  of  data  at  WMQ  at  a  local  distance  frcxn  foe  epicenter  illus¬ 
trate  foe  use  of  beta-stations,  even  if  foe  nature  of  foe  source  in  this  case  perhaps  could  be 
established  with  sufficient  confidence  from  only  teleseismic  data. 

The  value  of  identification  of  depth  phases  for  constraining  focal  depth  is  well 
known  and  demonstrated  in  this  case.  The  capability  of  arrays  to  estimate  direction  of 
approach  of  prinuuy  as  well  as  secondary  phases  enhances  foe  confidence  in  foe  identifica¬ 
tion  of  secondary  phases  as  depth  phases. 
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COMPARISON  OF  SHORT  PERIOD  WAVE  FORMS 


20  30  40  50 

TIME  (s) 

Figure  16:  Comparison  of  waveforms  at  NRAO  of  events  at  and  near  the  Lop  Not  test 

base  recorded  at  the  IMS  array  NRAO. 
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The  earthquakes  with  q>icenier  near  the  Lop  Nur  test  base  enyhasizes  the  value  of 
discrimination  case  studies  for  this  area  with  closely  spaced  e]q)losions  and  earthquakes. 
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Improvement  of  Regional  Seismic  Event 
Locations  for  Global  Monitoring 

Cari  R(»nney,  Vladislav  Ryalx^ 


Summary 

This  paper  investigates  improvements  in  regional  travel-time  curves,  and  die  identification 
of  regicMial  phases,  that  can  be  derived  from  Imig-range  seismic  refraction  profiles,  togetiier  with 
known  mine  blasts.  Deqi  Seismic  Sounding  (DSS)  observations  within  the  Baltic  ^eld  were 
used  to  derive  a  new  mean  regional  travel-time  model  for  Pn-waves  tiiat  extoids  to  distances  of 
1,400  km.  However,  pronounced  variations  of  Pa-wave  prc^gation  parameters  were  observed 
within  this  region,  iqiparently  caused  by  lateral  structural  inhomogeneities  of  the  earth's  crust  and 
upper  mantle.  Fw  example,  Pn  arrival  times  in  southern  Nmway  and  Sweden  at  distances  from 
200  to  SOO  km  are  approximately  2-3  sec  earlier  than  in  southern  Finland.  Consequently,  nuq;>s  of 
Pn  local  travel-time  conections  with  reflect  to  the  new  mean  regional  Pn  curve  were  cmxqnled 
for  the  NORESS,  ARCESS,  and  FINESA  arrays.  Comparison  of  the  Intelligent  Monitoring  Sys¬ 
tem  (IMS)  locations,  using  the  standard  IMS  r^erence  travel-time  curves,  with  those  d^ermined 
using  the  new  mean  regional  travel-time  curves  together  with  the  local  corrections,  showed  tiiat 
relocated  events  could  differ  with  the  IMS  locations  by  20  km,  and  the  difference  in  origin  time 
was  as  large  as  ±1.5  sec.  Simplifying  assumptions  in  the  IMS  about  the  travel-times  and  velocity 
of  Lg-waves  at  near-regional  distances  have  led  to  distance-dqiendent  location  biases  20  to  30 
km  for  mine  blasts  at  known  locations.  Using  DSS  data  to  help  define  the  pattern  tiie  observed 
waves  in  the  northern  Baltic  region,  we  found  that  the  assumed  Lg  waves  were  actually  SmS 
(shear-wave  reflection  from  the  Moho),  and  applying  the  appropriate  travel-time  curve  eliminates 
the  distance  component  of  this  bias.  Our  analysis  showed  that  large  travel-time  variatimis,  of  tiie 
order  of  10  sec  for  Pn-waves,  are  observed  in  different  regions  of  the  world.  We  came  to  the  con¬ 
clusion  that  none  of  the  existing  standard  global  travel-time  curves  can  cmiristently  produce 
regional  event  location  accuracies  adequate  to  meet  global  nuclear  test  ban  monitoring  objectives. 
However,  a  composite  Pn  travel-time  curve  that  we  define  for  Precambrian  shields  appears  to  be 
transpOTtable  to  other  similar  structural  provinces  of  tiie  worid,  which  corr^irise  much  of  the  con¬ 
tinental  area.  We  propose  that  this  curve  be  applied  in  such  provinces  for  monittxing  purposes. 
Tectonically  active  regions  will  probably  require  individual,  region-by-regimi  calibration. 
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Introduction 


The  effectiveness  of  a  seismic  network  in  a  cominehensive  nuclear  test  ban  moaitoring 
role  may  be  critically  dependent  on  the  accuracy  with  which  it  is  able  to  locate  seisnac  events. 
Get^hysical  and  other  techniques  designed  to  pinpoint  the  cavity  from  a  concealed  underground 
explosion,  <x  to  obtain  jnoof  of  its  nuclear  ori^n,  have  been  found  to  be  infective  when  qiplied 
to  searching  regions  of  the  size  of  normal  seismic  epicentral  uncertainty  areas  (Romney,  1990). 
F(»r  seismic  locations  using  just  a  few  stations,  as  will  be  tire  case  for  numerous  low  magnitude 
ev«its  near  the  network’s  detection  threshokl,  hi^y  accurate  regional  travel-time  tables,  phase 
timing,  and  azimuth  estimation,  as  well  as  the  ctxrect  identification  of  phases,  are  all  crudal  in 
iq)proaching  the  desired  accuracies.  The  objective  of  this  paper  is  to  show  tiutt  inqnoved  location 
capability  can  result  from  the  application  of  pn^gation  parameters  of  major  regional  phases 
derived  mainly  frmn  deep  seismic  sounding  (DSS)  profiles  within  the  Baltic  shield  and  neighbor¬ 
ing  areas.  We  will  use  ^ta  from  the  Intelligent  Mcmitt^g  System  (IMS)  [Bache  et  aL,  1S190] 
and  its  associated  regional  seismic  stations  in  Scandinavia  to  model  a  moiutoring  network  and  to 
denoonstrate  that  improvements  will  result  from  incotpmating  region-specific  propagation  param¬ 
eters. 


Event  locations  in  general  have  two  sources  of  error,  random  and  systematic.  One  souroe 
of  random  errors  is  measurement  error,  which  is  usually  a  function  of  signal-to-noise  ratio(SNR). 
Systematic  errors,  however,  are  caused  by  misidentification  of  phases,  inaccurate  travel-time 
titles,  and  biased  azimuth  estimates  [Bratt  and  Bache,  1988;  Myklwltveit  et  al^  1990].  The  most 
direct  application  of  DSS  observations  and  mine  blasts  with  well  determined  locati(»is  as  a  ttxd 
for  improving  event  locations  is  to  reduce  the  systematic  errors  due  to  inaccurate  travel-time 
tables,  as  well  as  to  the  misidentification  of  regional  phases.  Both  of  these  triplications  will  be 
investigated  in  diis  paper. 

Travel-time  tables  used  for  locating  local  and  regional  seismic  events  ate  generally  based 
on  one-dimensional  velocity  models.  For  the  IMS,  this  model  consists  of  two  crustal  and  two 
upper  mantle  plane  homogeneous  layers  [Bache  et  al.,  1990;  Bratt  et  al.,  1990].  The  crustal  thick¬ 
ness  for  this  model  is  40  km,  and  the  Moho  depth  and  the  layer  velocities  ate  an  average  fit  to  the 
region  under  consideration.  There  are,  however,  large  deviations  from  this  model  within  die  Bal¬ 
tic  shield.  According  to  detailed  DSS  observations,  the  Moho  dqith  varies  firom  approximately  30 
km  in  southern  Norway  and  southern  Sweden  to  about  SS  km  and  greater  in  southern  Hnland  and 
some  other  areas  [Luosto,  1991;  Sharov,  1991;  Ryaboy,  1990;  Kinck  et  aL,  1990;  BABEL,  1993a, 
1993b].  The  crustal  thickening  is  mainly  due  to  a  layer  with  a  F-wave  velocity  of  7.0-7.4  km/sec 
in  the  lower  crust  [Koriionen  et  al.,  1990;  Kinck  et  al.,  1990].  Travel-time  tables  used  fm:  locating 
events  usually  do  not  take  into  consideration  such  regionally  variable  velocity  structures.  The 
result  is  location  bias  in  areas  where  the  structure  departs  from  the  regional  mean. 

Misidentification  of  phases  can  have  an  even  larger  effect  Variations  of  tegioiuil  phase 
travel-times  and  amplitudes  make  phase  identification  a  difficult  task,  especially  where  informa¬ 
tion  on  apparent  surface  velocity  is  limited  or  subject  to  uncertainty  as  in  the  case  of  small  arrays, 
or  not  available  at  all  as  in  the  case  of  single  stations.  Because  DSS  data  are  observed  from  con¬ 
trolled  sources  and  sampled  at  numerous  fairly  regular  distance  intervals,  and  on  reversed  and 
overlapping  profiles,  the  entire  wave  pattern  of  initial  and  secondary  arrivals  becomes  visible  and 
regional  phases  can  be  identified  more  reliably. 
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In  what  follows,  data  from  DSS  profiles  in  the  Baltic  shield  are  applied  to  these  two  prob¬ 
lems:  namely,  defining  regionally  and  locally  corrected  travel  time  curves,  and  inqnoving  phase 
identification.  While  this  weak  is  primarily  focused  on  Scandinavia  and  the  study  of  IMS  data, 
the  methodology  can  also  be  applied  to  data  from  other  regional  seismic  networks.  In  later  sec¬ 
tions  of  this  paper  we  describe  some  ideas  on  extending  diis  methodology  to  global  monitoring. 

DSS  Observations 

Scandinavia  and  northwestern  Russia  have  been  covered  with  a  dense  network  of  DSS 
profile  observations  crossing  this  region  in  different  directions  (Figure  17).  The  results  of  these 
studies  have  been  published  in  numerous  papers,  monographs,  and  dissertations  [e.g.  Proceed¬ 
ings...,  1971,  1990;  EUGENO-S...,  1988;  Guggisberg,  1986;  Stangl,  1990;  Study  of  die  deq> 
structure...,  1986;  Ryaboy,  1989;  Structure...,  1991;  Imaging  and...,  1991;  Luosto  et  al.,  1984, 
1989, 1990;  Korbonen  et  al.,  1990].  The  lengths  of  the  DSS  profiles  vary  finom  200  to  2000  km, 
and  they  provide  observations  of  different  kinds  of  P-  and  S-  waves  travelling  in  the  earth’s  crust 
and  upper  mantle.  Field  observations  were  usually  carried  out  along  reversed  and  overlapping 
lines  for  dip  control  and  continuity.  Seismic  signals  were  typically  recorded  in  the  fi«quency 
band  from  1-2  Hz  up  to  20-25  Hz  by  vertical  and,  recently,  horizontal  seismometers. 

We  collected  and  analyzed  the  available  DSS  P-wave  record  sections  and  published 
travel-time  observations,  as  well  as  data  from  recordings  of  chemical  explosions  with  weU  known 
locations  and  origin  times  for  the  Baltic  shield  [Dahlman,  1967;  Porkka  et  al.,  1972;  Sellevoll  and 
Pomeroy,  1968;  Bath,  1981;  Cassell  et  al.,  1983].  A  number  of  shots  detonated  during  the  1979 
FENNOLORA  DSS  projec.t  from  nine  major  shot  points  were  recorded  by  the  NORSAR  array  at 
regional  distances  [Mereu  et  al.,  1983].  Record  sections  were  constructed  for  all  of  the  ^ot 
points  using  filtered  beams  at  each  of  the  seven  NORSAR  sub-arrays.  The  higher  coherency 
within  the  subarrays  improved  the  SNR  for  these  beams  relative  to  beams  constructed  as  stacks  of 
all  NORSAR  sensors  [Ryaboy,  1990]. 

NORESS  recordings  of  shots  detonated  along  the  POLAR  DSS  profile  in  northern  Nor¬ 
way  and  Finland  south  of  ARCESS  [Luosto  et  al.,  1989]  were  also  anriyzed.  They  had  poor 
SNR,  and  beams  constructed  for  different  explosions  did  not  show  marked  improvement  in  SNR. 
The  SNR  was  improved  to  the  point  where  the  Pit-wave  could  be  identified  only  after  stacking 
beams  for  several  explosions  detonated  at  the  same  shot  point. 

Pit  waveforms  observed  in  this  region  usually  have  a  complicated  structure  at  distances 
greater  than  300  km.  This  is  observed  most  clearly  when  conducting  detailed  observations  along 
reversed  and  overlapping  lines.  The  emergence  of  Pit  as  the  first  arrival  is  marked  by  an  increase 
in  apparent  velocity  from  6.0-7.0  km/sec  (for  the  Pg  and  Pb  waves)  to  7.5-8.0  km/sec  or  higher. 
Travel-time  curves  of  Pit  first  arrivals  are  not  linear,  and  there  is  a  tendency  toward  a  rise  in  the 
apparent  velocities  with  increasing  distance  to  8.5-9.0  km/sec  at  ranges  of  1,(X)0  to  1,500  km. 
Strong  arrivals  with  high  apparent  surface  velocities  are  often  detected  as  later  phases  in  the  P 
wave  group.  The  interval  over  which  they  can  be  correlated  fluctuates  from  tens  to  hundreds  of 
km.  Large  later  arrivals,  as  a  rule,  are  characterized  by  lower  frequencies  and  at  distances  less 
than  the  PgIPn  crossover  point,  by  higher  apparent  surface  velocities  than  the  low-amplitude 
waves  recorded  as  first  arrivals. 
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In  the  record  sections  constructed  for  the  FENNOLORA  and  OHL  profiles  [Guggisberg, 
1986;  Ryaboy,  1990;  Hauser  and  Stangl,  1990],  a  shift  of  1.5  -  2.5  sec  was  observed  in  the  travel¬ 
time  curve  between  800-900  km,  with  the  first  arrival  decaying  rapidly  and  a  later  phase  appear¬ 
ing  as  the  first  arrival.  Beyond  900  km,  the  apparent  surface  velocity  of  the  first  arrivals  sharply 
increased  to  8.6-8.8  km/sec  accompanied  by  a  decrease  in  frequency  (Figure  18).  This  feature  of 
the  Pn  record  sections  is  typical  for  the  Baltic  shield  and  is  most  likely  caused  by  a  low  veloci^ 
layer  in  the  upper  mantle  at  depths  greater  than  1(X)  km  [Hauser  and  Stangl,  1990;  Ryaboy,  1990]. 
Attenuation  of  Pn  first  arrivals  and  smaller  offsets  were  also  observed  at  distances  less  than  800 
km  for  several  DSS  profiles  within  Scandinavia  [Lund,  1979;  Proceedings...,  1990].  Local  small- 
scale  inhomogeneities  and  thin  low-velocity  layers  in  the  upper  mantle  are  the  suspected  cause  of 
these  variations. 

Pn  propagation  parameters  such  as  travel-time,  cross-over  distance  with  Pg,  amplitudes, 
firequencies,  and  polarization  characteristics  are  variable  within  this  region.  For  example,  we 
compare  two  DSS  profiles:  a  North  Sea  -  southern  Norway  profile  [C^assell  et  al.,  1983]  and  the 
SVEKA  profile  [Luosto  et  ai,  1984].  The  latter  was  observed  in  southern  areas  of  Finland  and 
the  former  in  southern  Norway  from  land  recordings  of  offshore  explosions.  The  waveforms 
beyond  450  km  for  the  northern  end  of  the  North  Sea  -  southern  Norway  profile  were  recorded  at 
the  NORSAR  array. 

There  are  important  differences  between  seismic  wave  propagation  parameters  for  these 
two  profiles  (Figure  19).  Pn  travel-times  in  southern  Norway  are,  for  the  same  distances,  inoie 
than  3  sec  early  in  comparison  with  those  in  southern  Finland.  The  Pg  /  Pn  cross-over  point  is 
located  at  a  distance  of  150  km  in  southern  Norway,  and  270  km  in  southern  Finland.  The  first 
arrivals  on  the  record  section  for  the  SVEKA  profile  from  200  to  270  km  are  lower  crustal 
refracted  waves  with  apparent  velocities  of  7.35  lan/sec.  Arrivals  with  similar  apparent  velocities 
were  not  recorded  in  southern  Norway.  The  map  in  Rgure  17  includes  an  overlay  of  Pg  /  Pn 
cross-over  distances  constructed  from  DSS  profile  c  bservations.  This  distance  varies  within  the 
Baltic  shield  from  130  to  270  km. 

Defining  the  Mean  Regional  Scandinavian  (MRS)  Pn  IVavel-Time  Curve 

The  available  P/z-wave  record  sections  and  travel-time  curves  constructed  for  DSS  profiles 
described  above  were  collected  to  correlate  and  measure  travel-times  of  Pn  first  arrivals.  Obser¬ 
vations  for  the  FENNOLORA  DSS  profile  extended  to  a  distance  of  2000  km,  but  first  arrivals 
could  be  reliably  picked  only  up  to  distances  of  1300-1400  km  because  of  poor  SNR. 

Travel-time  measurements  were  made  from  the  Pn  record  sections  (measured  at  10  km 
spacing)  and  published  Pn  travel-times  were  also  used.  These  data  were  applied  to  produce  the 
composite  Pn  travel-time  plot  shown  in  Figure  20.  This  figure  is  based  on  analysis  of  approxi¬ 
mately  10,000  km  of  DSS  profiles  and  includes  1420  measurements.  The  observations  have  a  scat¬ 
ter  ranging  up  to  about  three  seconds  for  a  fixed  distance.  This  scatter  is  much  larger  than  the 
observational  error  and  thus,  is  likely  caused  by  factors  such  as  lateral  structural  variations  of  the 
earth’s  crust  and  upper  mantle.  Also  evident  in  the  data  is  a  discontinuity  near  850  km,  indicative 
of  a  low-velocity  zone. 
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Figure  18:  Pn-wavc  record  sections  for  northern  (top)  and  southern  (bottom)  areas  of  the  Bal¬ 
tic  shield.  The  northern  record  section  was  constructed  for  shot  point  I  of  the 
FENNOLORA  DSS  profile  (adapted  from  Hauser  and  Stangl,  1990).  The  southern 
record  section  is  from  the  Oslo-Helsinki-Leningrad  (OHL)  profile  overlain  by  the¬ 
oretical  travel-time  curves  of  waves  refracted  in  two  upper  mantle  layers  separated 
by  a  low  velocity  zone  (adapted  from  Ryaboy,  1990).  Note  the  delayed  travel-time 
branch  beyond  850-900  km,  followed  by  an  increased  apparent  velocity. 
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Comparison  of  P-wavc  record  sections  from  the  Nonh  Sea-southera  Norway  profile  (top)  with  the  SVEKA 
profile  in  southern  Finland  (bottom).  DSS  profiles  are  adapted  from  Cassell  et  aL,  (1983);  and  Luosto  et  al., 
(1984),  respectively.  Note  the  significant  differences  in  the  reduced  travel-times  of  Pn  (about  5.5  -  6.0  sec  in 
southern  Norway  and  about  8.5  -  9.0  sec  in  soudiem  Finland),  and  in  the  distance  to  the  Pg/Pn  cross-over 
points  (about  150  km  in  souths  Norway,  and  270  km  in  southern  Hnland). 
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plots  of  Pn  travel-times  (dots)  shown  in  Figure  20  (top).  Boxes  were  calculated  for  100  km  distance  windows.  The  hor¬ 
izontal  line  in  the  interior  of  each  box  is  located  at  the  median  of  the  data.  The  top  and  bottom  of  each  box  represent  the 
second  and  third  quartiles  of  the  data.  The  dotted  vertical  lines  extend  through  the  6rst  and  fourth  quartiles  to  the 
extreme  values  of  the  data.  The  solid  line  labeled  MRS  is  the  proposed  new  mean  regional  Scandinavian  reference  Pn 
travel-time  curve. 


Three  travel-time  curves  are  superin^msed  on  the  observed  DSS  travel-time  data  shown  in 
the  upper  section  of  Hgure  20.  Two  of  these,  IMS  and  HEL,  are  the  curves  used  to  locate  events 
for  the  IMS  and  Helsinki  bulletins,  respectively.  Both  of  these  Pn  travel-time  curves  are  faiily 
good  fits  to  the  DSS  data  for  distances  less  than  SOO-600  km.  However,  beyond  about  600  km, 
DSS  observations  axe  generally  early  with  respect  to  the  IMS  Pn  reference  travel-time  curve  for 
the  Baltic  shield.  On  the  other  hand,  DSS  observations  at  these  longer  ranges  are  generally  late 
with  respect  to  the  Pn  travel-time  curve  currently  in  use  to  locate  events  for  the  Ifelsinid  bulletin 
(HEL).  The  third  travel-time  curve  on  Hgure  20,  labeled  MRS,  represents  our  current  best  esti¬ 
mate  of  a  generalized  mean  regional  Scandinavian  Pn  travel  time  curve  for  die  entire  region 
[Ryaboy,  1993].  Local  corrections  to  this  curve  will  be  discussed  subsequently. 

Our  preferred  model  to  fit  the  observed  DSS  Pn  travel-times  is  based  on  a  one-dinunsional 
spherical  earth  model  with  two  crustal  layers  and  a  low  velocity  zone  in  the  upper  mantle  (Hgure 
21).  The  MRS  travel  time  curve  shown  on  Hgure  4  is  calculated  from  this  model.  Wocity  noodels 
underlying  the  IMS  and  HEL  travel-time  curves  are  also  shown  on  Figure  21.  Our  preferred  model 
is,  of  course,  not  unique,  and  is  only  one  representative  of  a  family  of  velocity  models  consistent 
with  the  observed  Pn  travel-times.  However,  it  does  have  the  merit  of  providing  a  better  fit  to  die 
DSS  Pn  observations  than  do  the  other  two  models,  as  shown  by  its  excellent  fit  to  the  box  plots 
on  the  lower  section  of  Figure  20,  representing  the  DSS  observations  grouped  within  1(X)  km  dis¬ 
tance  windows. 

Sub-Regional  and  Local  TVavel-Time  Corrections 

We  interpret  the  observational  deviations  from  the  MRS  curve  (Hgure  20),  which  have  a 
standard  deviation  of  ±0.65  sec,  to  be  mainly  the  result  of  sub-regional  structural  and  velocity  inh¬ 
omogeneities.  It  was  earlier  noted  that  the  PgIPn  crossover  distance  varied  systematically  over  this 
region.  This  observation  is,  of  course,  a  reflection  of  systematic  travel  time  variations.  In  southern 
Norway  and  southern  Sweden  the  F/i-wave  is  recorded  1.0-1.5  sec  early,  and  in  southern  Hnland 
1.0- 1.5  sec  late,  relative  to  our  preferred  Pn  travel-time  curve.  The  earliest  Pn  travel-times  were 
observed  within  a  narrow  zone  of  the  Oslo  graben  (50-100  km  wide)  south  of  the  NORSAR/NOR- 
ESS  arrays.  In  the  southern  part  of  the  FENNOLORA  profile,  at  distances  of  800-1000  km,  the 
P/i-wave  is  recorded  1.5-2.0  sec  earlier  than  in  the  northern  part  of  the  profile.  This  difference  in 
Pn  travel-times  cannot  be  explained  by  variations  of  crustal  structure  along  the  profile,  and  proba¬ 
bly  indicates  a  zone  of  increased  velocity  in  the  lower  lithosphere  beneath  southern  Sweden.  This 
zone  was  detected  as  a  result  of  a  joint  interpretation  of  seismic  and  gravity  observations  for  the 
Oslo-Helsinki-Leningrad  profile  [Ryaboy,  1990],  and  is  supported  by  observations  of  teleseismic 
data  [Husebye  et  al.,  1986].  The  fact  that  these  variations  appear  to  be  systematic  over  subregions 
of  significant  size  suggests  that  it  is  feasible  to  take  them  into  account. 

The  method  to  be  investigated  here  for  estimating  local  travel-time  corrections  for  improv¬ 
ing  automatic  event  locations  takes  advantage  of  the  extensive  work  throughout  the  entire  region 
of  interest  for  determining  the  depth  of  the  Moho.  The  residuals  of  the  observed  Pn  travel-times 
with  respect  to  the  MRS  curve  (TT,qs)  are  highly  correlated  with  variations  of  the  Moho  depth  (Fig¬ 
ure  22)  (coefficient  of  linear  correlation  of  0.86  and  standard  deviation  of  ±0.39  sec): 
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Figure  21:  Reference  velocity  models  of  the  earth  crust  and  upper  mantle  for  the  Baltic  shield.  Numbers  are  P-wave  velodties  in 
km/sec,  with  layer  thidoiess  (km)  in  parentheses.  The  model  on  the  right,  based  on  DSS  data,  results  in  die  h^tS 
travel-time  curve. 
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Figure  22;  Pn  travel-time  residuals  (with  respect  to  the  MRS  curve)  vs.  the  Moho  depth  tot  the  Baldc  shield.  The  residuals  are 
highly  correlated  with  variations  of  the  Moho  d^th.  This  empirical  relation  and  a  Moho  depdi  map  were  used  to  calcu¬ 
late  ray-path  dependent  local  corrections  of  Pn  ttavd-dmes  for  regional  event  locations  within  the  Baltic  shield  and 
neighboring  areas. 


TTre,  -  -4.802  +  0.058  •  (M^  +  Mgec). 


(1) 


wbae  Mq,  and  M,ec  are  the  Moho  dq)ths  in  km  beneath  shot  pmnt  and  receiver  on  the  DSS  pro¬ 
files.  This  is  an  expected  result,  ofcourse,  since  these  identical  DSS  trevel-diDes  entered  into  the 
estimalesof  MdiodepdL  Values  were  measured  fiom  die  smoodiedTTKs  vs.  distance 

curves  constructed  for  different  shot  ptunts  on  DSS  profiles. 

The  empirical  relatioo  (1)  and  die  Moho  dqith  map  were  used  to  calculate,  on  a  r^ularly 
spaced  grid,  the  ray-path  dqiendent  local  conections  mi^y  caused  by  lateral  variations  of  the 
Mtdio  depth.  A  Medio  depth  nuqi  for  the  Baltic  shield  and  net^boring  r^ions  has  been  con¬ 
structed  based  on  the  new  and  most  reliable  published  data  [BABEL,  1993a,  1993b;  IQnck  et  aL, 
1990;  Kmhonen  et  al.,  1990;  Luosto,  1991;  Sharov,  1991;  Structure  and  dynamics...,  1991; 
of  the  deep...,  1986].  The  crustal  thickness  varies  from  approximately  30  km  in  soudiem  Norway 
to  50-55  km  and  more  in  soudiemHnland  and  Sweden.  The  Moho  map  was  digitized  and  included 
in  a  gee^hysical  database  at  Cornell  University  [Helding  era/.,  1993].  This  was  then  used  to  cal¬ 
culate  the  corrections  to  the  mean  regional  travel-time  curve  throughout  die  region  and  for  eadi 
station.  As  an  exanqiie,  Hgure  23  shows  the  noap  of  local  correctitxis  for  die  FINESA  station.  The 
maps  for  NORESS  and  ARCESS  are  timilai,  with  differences  caused  by  the  crustal  thickness  vari¬ 
ations  beneath  these  stations.  Pn  local  cmrections  for  Scandinavia  range  within  approxunaiely 
±1.5  sec. 

However,  an  improvement  in  the  Pn  tables  alone  is  not  enough  to  give  dramatic  results  for 
distances  less  than  800-1000  km  within  the  Baltic  shield.  Based  on  our  experience  widi  automatic 
^tems  for  regional  event  location,  we  have  fmind  that  event  mislocation,  in  some  distance 
ranges,  is  as  much  a  function  of  assumptions  about  regional  phase  identification  as  it  is  (tf  r^ional 
errors  in  travel-time  tables.  In  the  following  section  of  this  pr^ier,  the  problem  of  phase  identifica¬ 
tions  fOT  regional  event  locations  is  investigated. 

Improved  Phase  Identification 

A  detailed  study  of  ARCESS  epicenters  in  the  nexthem  Baltic  shield  from  mine  blasts  hav¬ 
ing  known  locations  demonstrated  bias  in  distance  and  in  azimuth  [Suteau-Henson,  IS^l;  Sutean- 
Henson  et  al.,  1992].  We  consider  in  this  paper  only  the  bias  in  distance.  To  study  tins  phenome¬ 
non,  we  carefully  analyzed  selected  ARCESS  reemdings  of  mine  blasts  from  four  known  locations: 
Varanger  (V),  Saattopora  (S),  Kiruna  (K),  and  Elijaervi  (E)  (Hgure  24).  Distances  to  the  approxi¬ 
mate  centers  of  these  mines  from  ARCESS  are  156  km,  200  km,  287  km,  and  421  km,  respectively. 
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Figure  24:  Location  map  of  the  ARCESS  array  (mangle)  and  seismic  events  (small  squares) 
from  mine  blasts  recorded  within  the  northern  Baltic  shield  and  listed  in  the  IMS 
bulletin.  Large  circles  show  distance  from  ARCESS  at  100  km  intervals.  The 
Varanger,  Saattopora,  Kiruna.^  and  Elijaervi  mines,  studied  in  this  w<»ic,  are  desig¬ 
nated  by  open  circles  and  capital  letters  V,  S,  K,  and  E,  respectively.  Saattopora 
mine  blast  locations  at  200  km  do  not  have  a  bias  in  distance.  Kiruna  and  Elijaervi 
mine  blast  locations  at  distances  greater  than  200  km  are  biased  toward  ARCESS, 
while  Varanger  locations  at  a  lesser  distance  are  biased  away  from  ARCESS. 
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Oon^Miing  the  mine  kx:adons  to  IMS  ^icenters,  adistanoe-dqiendent  trnid  in  the  svenife 
bias  was  noted,  as  shown  in  die  following  Ihble  S: 

TsUeS: 


Mine 

Mine 

Distance 

frxxn 

ARCESS. 
in  km 

Avenge 
Epicentral 
Distance, 
in  km 

Bias, 
in  km 

Varanger 

156 

171 

+  15 

Saattt^iora 

0 

Kiruna 

287 

272 

-15 

Elijaervi 

421 

391 

-30 

The  huger  distance  biases  are  much  greater  than  would  be  expected  finom  the  variations 
caused  by  local  or  sub-regional  perturbations  in  the  travel-time  curves,  especially  since  DSS  data 
in  the  regions  of  the  mines  show  that  P/i-wave  travel-time  variations  are  only  appmcimaiely  0  J 
second,  comparable  to  the  measurement  error  (Hgure  23).  Thus,  the  bias  is  lik^y  die  result  of  mis¬ 
timing  <x  misidentification  of  phases.  We  turn  to  DSS  observations  in  Scandinavia  to  help  resolve 
this  problem. 

Detailed  observations  based  on  recordings  of  P-  and  5-waves  were  earned  out  in  this  area 
alcHig  the  FENNOLORA  [Guggisbeig,  1986;  Hauser  and  Stangl,  1990;  Stangl,  1990]  and  POLAR 
[Luosto  et  al.j  1989]  long-range  DSS  profiles  (Figure  17).  As  an  exanqile,  Hgure  25  shows  die 
record  sections  of  P-  and  5-waves  for  the  POLAR  DSS  profile.  Pg-  and  5g-waves,  refracted  in  die 
upper  crustal  layer,  were  observed  as  first  arrivals  within  the  P-  and  5-wave  groups,  reflectively, 
up  to  a  distance  of  120-140  km.  Pb-  and  5b-varves,  refiacted  in  die  lower  cnist,  were  usually 
observed  as  first  arrivals  up  to  that  distance  and  continuing  to  about  200  km,  where  Ph  and  Sn 
emerge  as  first  arrivals.  High  angilitude  PmP  and  SmS  reflected  waves  finm  the  Moho  were 
observed  in  the  later  arrivals  of  the  P-  and  5-wave  groups,  reflectively,  at  distances  fiom  fipiox- 
imately  100  to  350-400  km.  \^thin  this  range  the  an^ilitodes  of  these  reflected  waves  netypicaUy 
much  larger  than  the  refracted  first  arrivals. 


The  observed  travel-times  of  Pg,  Pb,  Pn,  Sg,  Sb,  and  5/i-waves  recorded  as  first  arrivals 
along  the  POLAR  DSS  profile,  as  well  as  the  travel-times  of  the  most  intensive  waves  (PmP  and 
SmS)  recorded  as  late  arrivals,  are  shown  on  Hgure  26  and  Hgure  27.  Also  shown  on  these  fig¬ 
ures  are  the  IMS  reference  travel-time  curv^  for  the  first  arrivals,  along  with  curves  for  PmP  and 
5m5-waves  derived  from  the  DSS  observations.  The  observed  P-wave  arrival  times  fixxn  the 
POLAR  profile  fit,  in  general,  the  IMS  reference  travel-time  curves  fte'  Pg,  Pb  and  Pn-waves. 
However,  Pg-waves  were  not  recorded  along  the  POLAR  profile  at  distances  beyond  the  cross¬ 
over  point  (Figure  25).  Observed  arrival  times  of  Sg  and  5b  also  fit  the  IMS  travel-time  curve  for 
Lg.  Sg  was  not  recorded  on  the  POLAR  profile  at  distances  greater  than  about  200  km  nor  was 
Lg.  5ii-waves  recorded  along  the  POLAR  profile  fit  the  reference  IMS  curve  between  the  cross- 
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Figure  25:  Trace-normalized  record  sections  of  S-  (top)  and  P-waves  (bottom)  recorded  on  the  POLAR  DSS  profile  fixxn  shot 
point  A  (adapted  from  Luosto  era/.,  1989).  Band-pass  filters  of2-.6  Hz  and  2-15  Hz  were  lulled  to  die  reomds  of  5- 
and  P-waves,  re^tectively.  It  may  be  seen  that  the  most  energetic  arrivals  betweoi  about  80  km  and  280  km  are  PmP 
and  SmS. 


P-WAVE  TRAVEL-TIMES  FOR  POLAR  DSS  PROFILE 
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S-WAVE  travel-times  FOR  POLAR  DSS  PROFILE 
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over  point  and  a  distance  of  about  300  km.  Beyond  300  km,  Sn  was  observed  to  arrive  systemati¬ 
cally  earlier. 

We  note,  that  PmP  and  5m5-waves  were  not  included  in  the  initial  list  of  major  regional 
phases  for  the  IMS,  but  these  waves  can  be  easily  incorporated  and  used  for  event  locations. 
Indeed,  the  IMS  was  specifically  designed  to  incoiporate  new  knowledge  of  this  kind  as  it 
emerges  fiom  research.  Examples  of  using  PmP-wave  for  regional  event  locations  were  recently 
published  [Quin  and  Thurber,  1992]. 

Figure  28  and  Figure  29  show  waveforms  recorded  at  ARCESS  of  selected  mine  blasts 
from  the  Varanger  and  Kiruna  mines.  Phase  identifications  determined  by  the  IMS  using  an  auto¬ 
mated  procedure  and  after  review  by  an  analyst  are  marked  above  each  waveform,  ^^cal  lines 
running  through  all  traces  are  label^  at  the  top  with  the  phase  identification  based  on  DSS  data, 
using  travel  times  form  Figures  26  and  27. 

For  events  at  both  mines,  we  note  that  the  beginning  of  the  intensive  waves  within  the  S- 
wave  group  were  identified  in  the  IMS  mainly  as  Lg.  Despite  the  similarity  in  the  waveforms, 
there  is  a  considerable  amount  of  scatter  in  the  Lg  phase  piclu,  which  is  the  main  cause  of  the  dis¬ 
tance  scatter  in  the  event  locations.  Based  on  our  interpretation  using  DSS  data,  the  wave  labeled 
as  Lg  by  the  IMS  is  actually  SmS.  Difference  between  the  SmS  and  Lg  travel-time  curves  are  dis¬ 
tance  dependent  and  may  be  seen  to  explain  the  observed  distance  bias  for  event  locations  in 
Table  5.  For  example,  the  Lg  travel-time  curve  from  the  IMS  bulletin  is  approximately  3  seconds 
earlier  than  SmS  at  the  distance  of  the  Varanger  mine  blasts.  Interpreting  SmS  as  Lg,  therefore, 
results  in  a  biased  location  away  from  ARCESS.  At  the  distance  of  Scattopora,  about  200  km,  tire 
Lg  curve  passes  through  the  SmS  observations  and  thus  there  is  no  resulting  bias  from  interpreting 
SmS  as  Lg.  At  the  greater  distances  of  the  Kiruna  and  Elijaervi  mines,  the  Lg  curve  is  3  seconds 
and  (perhaps)  6  seconds,  respectively,  later  than  SmS,  resulting  in  biases  toward  the  station  if  SmS 
is  interpreted  to  be  Lg. 

Our  SmS  interpretation  of  the  wave  labeled  as  Lg  by  the  IMS  is  consistent  with  studies  of 
major  regional  phases  recorded  by  the  NORESS,  ARCESS,  and  FINES  A  arrays  within  the  Baltic 
shield  [VogQord  et  al.,  1991;  VogQord  and  Langston,  1990,  1992]  and  results  of  modelling  L^- 
wave  propagation  [Kennett,  1992a,  1992b].  Lg  is  believed  to  be  a  superposition  of  5-waves  post- 
criticdly  incident  on  the  Moho  and  multiply  reflected  in  the  crust  [Bouchon,  1982;  Campillo  and 
Paul,  1^2;  Kim,  1987].  The  Lg  group  velocity  is  near  3.S  km/sec  and  the  phase  has  predomiruuit 
frequencies  of  1-2  Hz  and  lower  [Press  and  Ewing,  1952;  Kennett  et  al.,  1985;  Ruzaildn  et  al., 
1977].  SmS  consists  of  pre-critical,  critical,  and  post-critical  reflected  shear  waves  from  the 
Moho.  Analysis  of  DSS  data  suggests  that  SmS  has  a  phase  velocity  larger  than  that  reported  for 
Lg  [for  a  review  of  Lg  phase  and  group  velocity  measurements,  see  e.g.  Blandford,  1981]. 

Figure  30  shows  the  apparent  velocity  for  the  phase  identified  as  Lg  by  the  IMS  as  a  func¬ 
tion  of  distance  from  the  ARC^S  array.  Also  shown  in  Figure  30  is  the  SmS  apparent  velocity 
curve  inferred  from  the  POLAR  DSS  profile.  The  SmS  curve  may  be  seen  to  be  a  good  fit  to  the 
IMS  data  within  the  distance  range  from  100  to  4(X)  km,  approximately.  At  distances  between  400 
and  LCXX)  km,  typical  Lg  apparent  velocities  (slightly  above  3.5  km/sec)  are  observed.  At  dis¬ 
tances  less  than  100  km,  however,  the  apparent  velocity  inferred  from  the  IMS  observations  is 
between  4.0-5.0  km/sec.  These  values  are  too  large  for  Lg  and  too  small  for  SmS,  but  are  typical 
for  reflected  waves  from  mid-crustal  boundaries  in  the  Bdtic  shield  [Luosto  et  al.,  1989;  Hauser 
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ARCESS  RECORDINGS  OF  VARANGER  MINE  BLASTS  (sz,  f-2-5  Hz,  array-mine  distance  Is  156  km) 

Pg  pnP  Sg(X.g)  SmS 
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channels).  Phase  identification,  indicated  above  each  trace  arc  based  on  the  IMS  bulletin.  Vertical  solid  lines  are  times 
of  arrivals  of  the  expected  major  P-  and  5-phases  {Pg,  PmP,  Sg,  and  5m5)  inferred  from  the  POLAR  profile  observations 
for  the  distance  of  156  km  between  ARCESS  and  the  Varanger  niine  (see  the  Figures  14  and  15).  Note  that  the  most 
prominent  P  arrivals,  PmP  was  not  used  in  the  IMS  analysis.  Sg  and  Lg  travel-time  curves  are  identical  at  this  distance. 
The  most  prominent  arrival  in  the  5  group  was  generally  interpreted  by  the  IMS  as  Lg. 


ARCBSS  RECORDINGS  OF  SELECTED  XIRUNA  MIME  BLASTS.  R-290  Xm 
(central  vertical  short  period  channels/  f-2-8  Hz) 
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channels).  Phase  identifications  indicated  above  each  trace  are  based  on  the  IMS  bulletin.  >^sitical  solid  lines  are  times 
of  arrivals  of  major  P-  and  S-phases  {Pn,  PmP,  Sn,  and  SmS)  inferred  from  the  POLAR  profile  observations  for  the  dis¬ 
tance  of  287  km  between  ARCESS  and  the  Kiruna  mine.  The  vertical  dashed  line  is  the  Lg-wave  arrival  time  accord¬ 
ing  to  the  IMS  reference  travel-time  tables.  PmP  is  not  detectable,  but  the  start  of  the  intensive  waves  in  the  5  group 
appears  to  be  SmS. 
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Figure  30:  Top:  Lg  apparent  velocity  vs.  distance  from  the  ARCESS  array  according  xof-k 
analysis  done  by  the  IMS  (dots).  Bottom:  box  plots  vs.  distance  for  the  sa^ 
IMS  data  (dots)  shown  above.  Boxes  were  calculated  for  100  km  distance  win¬ 
dows.  The  horizontal  line  in  the  interior  of  the  each  box  is  located  at  the  m^an 
of  the  data.  The  height  of  each  box  is  the  difference  between  the  third  quartile  of 
the  data  and  the  first  quartile.  The  dotted  vertical  lines  extend  to  the  extreme  val¬ 
ues  of  the  data.  Top  and  bottom:  the  solid  line  is  the  SmS  apparent  velocity  vs. 
distance  curve  inferred  from  the  POLAR  DSS  profile  observations;  the  dashed 
horizontal  line  indicates  3.5  km/sec,  the  Lg-wave  group  velocity. 
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and  Stangl,  1990].  These  arrivals  can  be  seen  on  the  DSS  profile  in  Hgure  25  (tc^)  at  distances 
less  than  100  km. 

Figure  31  shows  that  die  predominant  frequencies  of  the  phase  identified  asLg  by  the  IMS 
decreases  with  distance  from  8-10  Hz  at  100  km  to  2-4  Hz  at  400  km.  Beyond  400  km.  die  mean 
frequencies  approach  2  Hz  or  less,  which  is  more  typical  for  Lg  energy.  Since  die  lecocded  events 
with  only  a  few  exceptions  have  similar  magnitude,  and  one  can  see  a  sharp  frequency  decrease 
for  distances  greater  than  400  km,  variations  in  the  dominant  frequency  are  probably  not  due  to 
systematic  differences  in  the  source  size  or  attenuation.  High  frequencies,  such  as  those  observed 
at  distances  less  than  400  km,  are  atypical  for  Ig-waves  and  support  a  SmS  interpretation. 

We  conclude  that  the  first  high  amplitude  phases  recorded  .within  die  5-wave  group 
between  100-400  km  in  the  northern  Baltic  shield  and  having  apparent  velocities  of  approxi¬ 
mately  5.0  km/sec  to  4.0  km/sec  are  5/R5-waves,  and  the  high  amplitude  phases  recorded  at  dis¬ 
tances  greater  than  400-500  km  are  Lg-waves  having  apparent  velocities  of  qiproximately  3.5 
km/sec.  We  have  shown  that  mis-identification  of  SmS  as  Lg  creates  a  distance  dependent  bias  in 
event  locations  at  distances  from  approximately  100  km  up  to  400-500  km.  The  prmiounced  bias 
was  not  observed  at  about  200  km  because  this  is  the  cross-over  distance  for  the  Lg  and  SmS 
travel-time  curves. 

There  are  also  other  sources  of  phase  mis-identifications.  For  example,  we  have  noted  that 
PmP  on  ARCESS  recordings  of  the  Elijaervi  mine  blasts  have  been  mis-identified  as  Pg.  It  was 
also  recently  shown  that  GERES  S  recordings  of  Pg  have  been  mis-identified  as  Pn  [Grant  et  aL. 
1993].  Long-range  profile  observations  might  be  extremely  helpful  in  such  cases  for  improve¬ 
ment  of  phase  identifications. 

Relocations  of  the  Kiruna  mine  blasts 

To  test  the  influence  of  different  factors  on  regional  event  locations,  we  relocated  two 
groups  of  Kimna  mine  blasts:  the  first  located  by  data  from  ARCESS  alone  (25  events),  and  the 
second  by  data  from  two  or  three  arrays,  including  ARCESS  (29  events).  For  these  tests  we  used 
the  LocSAT  code  (S.  Bratt  and  W.  Nagy,  1991),  the  same  algorithm  that  was  used  for  the  IMS 
event  locations.  However,  we  substituted  the  new  MRS  travel-time  curve  for  tite  standard  IMS 
Pn  curve,  and  we  included  SmS  travel  times  from  Figure  27.  The  phase  at  ARCESS  labeled  Lg  by 
the  IMS  was  renamed  SmS,  but  all  other  phase  designations  and  the  arrival  times  were  unchanged 
from  the  IMS  picks.  Results  are  given  in  Table  6,  which  shows  that  the  distance  bias  is  essentially 
eliminated  for  locations  based  on  ARCESS  alone,  and  somewhat  reduced  for  locations  employing 
data  from  two  or  three  stations.  Our  expectation  is  that  application  of  the  same  care  in  interpreta¬ 
tion  of  data  from  the  other  stations  that  we  have  applied  to  ARCESS  data  would  result  in  further 
improvements  to  the  two  or  three  stations  locations. 
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Table  6:  DISTANCE  FROM  ARCESS  TO  KIRUNA  MINE  SHOT  LOCATIONS 


mean, 
in  km 

median, 
in  km 

S*d«f 

in±km 

ARCESS 

locations 

IMS  bull. 

270.84 

271.83 

12.11 

new  travel¬ 
time  tables 
aadSmS 

285.92 

287.25 

15.29 

2  or  3 
array 
locations 

IMS  bull. 

273.65 

276.76 

6.04 

new  travel¬ 
time  tables 
and  SmS 

278.95 

278.30 

6.63 

Note:  tnie  distance’ from  ARCESS  to  the  Kiruna  mine  is  287.3  km. 

We  should  explain  here  that  w  c  :  .ave  used  the  distance  of  the  epicenters  firom  ARCESS  as 
our  measure  of  improvement,  rather  ’ban  die  distance  to  the  mine,  because  there  is  also  known  to 
be  an  azimuthal  bias  attributable  to  small-scale  structural  variations  in  the  crust  and  upper  mantle, 
which  produces  mislocations  comparable  to  the  distance  bias  (Suteau-Henson  et  al,  1992)  (note 
also  Figure  24)  and  which  tends  to  obscure  the  improvements  in  distance  alone  that  are  under 
investigadon  here.  However,  auxiliaiy  tests  that  we  conducted,  using  azimuthal  corrections  pro¬ 
vided  by  Suteau-Henson,  showed  that  the  combined  effect  of  including  azimuthal  corrections, 
improved  travel  time  curves  and  improved  phase  identification  to  ARCESS  data  reduced  die 
mean  location  error  from  about  31  km  to  about  10  km.  Figure  32  shows  these  revised  locations;  it 
may  be  seen  that  both  the  azimuthal  and  distance  bias,  which  were  cleariy  evident  cm  Hgure  24, 
have  been  largely  removed.  As  previously  noted,  there  was  considerable  scatter  in  the  arrival 
times  of  Lg  (SmS)  as  picked  by  the  IMS,  which  iuxx>unts  for  some  of  the  remaining  location  scat¬ 
ter.  Figure  33  indicates  the  relative  effects  of  the  different  changes  made  in  the  tests  just 
described.  It  shows  that  the  most  significant  change,  by  far,  results  from  reinterpretation  of  Lg  as 
SmS;  by  itself,  this  eliminates  most  of  the  distance  bias  from  the  ARCESS  data.  The  new  Pn 
travel  time  data  further  reduces  the  bias,  and  an  independent  pick  of  the  arrival  times  of  the  strong 
shear  waves  --  with  the  knowledge  that  it  was  SmS  rather  than  Lg  -  gave  an  additional  improve¬ 
ment. 
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Figure  32: 


Map  of  Kiruna  mine  blasts  after  relocation.  Crosses  are  events  listed  in  the  IMS 
bulletin  and  relocated  using  the  new  travel-dme  tables  for  the  Baltic  shield  (for  the 
Pn-,  Sn-  and  5/n5-waves)  and  corrected  for  azimuthal  bias.  Large  circles  are  dis¬ 
tances  from  the  Kiruna  mine,  and  a  solid  line  shows  the  direction  from  the  Kiruna 
mine  to  ARCESS. 
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Figure  33:  Density  plots  of  distances  from  ARCESS  for  Kinina  mine  blast  locations  and  relo¬ 

cations.  Interpreting  the  strong  shear  waves  as  Lg  results  in  the  bias  evident  in  the 
"IMS  Bulletin"  plots;  re-interpreting  these  waves  as  SmS  corrects  most  of  this 
bias,  and  application  of  the  new  Pn  curves  and  more  consistent  measurements  of 
the  arrival  times  further  reduces  the  bias. 
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Comparison  of  regional  Pn  travel-time  curves  for  different  geological 
provinces 

Global  monitoring  for  nuclear  proliferadoa  is  expected  to  be  based  on  an  automated  sys¬ 
tem  that  includes  a  netwotk  of  high-quality  sdsmologi^  stations,  many  equqiped  with  anays, 
established  at  sites  having  substantially  different  cnistal  and  upper  mantle  structures.  To  assess 
the  influence  of  regional  structural  variations  on  travel  times  we  couqMted  the  new  MRS  Pn 
travel-time  curve  with  similar  curves  constructed  for  other  regions.  Hgm  34  shows  a  oouqiari- 
son  of  the  standard  global  Pn  travel-time  corves  titan  die  Jefitieys-BuUen,  Herrin  et  al^  and 
lASPEI-91  Seismological  Ihbles  [Kennett,  Engdahl,  1991]  with  several  regional  Pn  travel-time 
curves  inferred  tiom  explosion  seismology  observations.  Thtvel-time  curves  for  the  western 
USA,  Neath  American  platform,  and  Canadian  shield  were  constructed  by  smoothing  observa¬ 
tions  titan  nuclear  explosions  [Romney  et  at.»  1962]  and  chemical  explosions  in  Lake  Siqieiiar 
and  Hudstai  Bay  [Barr,  1967;  Hobson,  1967;  Mereu  and  Hunter,  1969;  Green  and  Hales,  1968; 
Masse,  1973].  The  Pn  travel-time  curves  for  the  Baltic  shield.  East  European  and  Siberian  plat¬ 
forms,  and  tia*  tectonically  active  regions  of  Hoi  Shan  and  the  Baikal  rift  zone  were  inferred  tiom 
DSS  profile  observations  and  tiom  nuclear  and  large  chemical  explosion  rectadings  by  seisnao- 
logical  networks  [\^imik,  Ryaboy,  1981;  Ryaboy,  1989;  Egorldn  et  al.,  1987]. 

The  curves  on  Figure  34  reveal  major  differences  between  travel-times  of  Pit-waves  prop¬ 
agating  in  tectonically  active  regions  and  travel-times  where  propagation  was  within  Precamhrian 
platforms.  It  may  be  noted  that  the  Pit  travel-time  curve  for  a  tectonically  active  region  of  the 
western  USA  is  characterized  by  substantially  greater  travel-times  (up  to  about  10  sec  within  die 
distance  range  tiom  S(X)  to  2000  km)  than  observed  within  Precambrian  geological  provinces.  Pit 
travel-times  for  intercontinental  tectonically  active  regions  (the  Hen  Shan  and  Baikal  rift  zone) 
are  approximately  5  sec  faster  than  fen*  the  western  USA  but  are  still  slower  than  Pit  travel-time 
curves  for  the  Precambrian  platforms  and  shields. 

These  Pit  travel-time  variations  provide  evidence  for  regions  of  large  scale  lateral  velocity 
inhtHnogeneities  in  the  earth’s  crust  and  upper  mande  with  linear  dimensions  of  lyiproximately 
1000-2(XX)  km.  The  regional  travel-time  curves,  of  course,  fit  different  geologic  provinces  only 
on  the  average.  To  use  these  curves  for  event  locations,  they  will  need  to  be  adjusted  to  specific 
areas  and  seismic  stations.  This  adjustment  means  that  path-dependent  corrections  sintilar  to 
those  in  Figure  23  will  be  needed,  based  on  available  geological  arid  geophysical  data,  to  account 
for  such  smaller  scale  features  as  deep  sedimentary  basins  and  changes  in  the  Medio  depth. 

We  note  that  the  standard  global  refnene^  Pn  curves  ~  the  lASPEI  1991  Seismological 
Tables,  the  Jeftieys-Bullen  Tables  and  the  Herrin  Tables  --  do  not  provide  a  good  fit  to  any  of  the 
observed  regional  curves  especially  at  distances  greater  than  500-600  km.  We  infer  that  none  of 
these  standard  global  reference  curves  will  produce  epicentral  location  accuracies  tixr  regional 
events  that  are  good  enough  for  global  monitoring  aims. 

Unfortunately,  accurate  regional  travel-time  curves  specific  to  many  of  the  regions 
important  for  global  monitoring  do  not  exist.  A  key  question  that  we  address  here  is:  what  would 
be  the  initial  assumptions  about  travel-times  for  such  regions?  Would  we  simply  apply  standard 
reference  curves  like  the  lASPEI  1991?  Or  can  we  do  better?  Fortunately,  there  is  much  infemna- 
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Figure  34;  Comparison  of  regional  Pn  travel-time  curves  for  different  geological  provinces  with  the  reference  Pn  travel-time 
curves  from  the  Jeffreys  and  Bullen,  Herrin  et  aL,  and  lASPEI  1991  Scismological  Tkbles.  Ancient  platforms  and 
shields  have  the  fastest  Pn  travel-times,  and  tectonically  active  regions  are  substantially  slower  (up  to  10  sec). 


tion  to  indicate  that  we  can  do  better,  at  least  under  some  circumstances. 
Transportability  of  a  Pn  TTavel-Tinie  Curve  for  Precambrian  Shields 


The  suiprisingly  good  agreement  among  the  Pn  curves  for  Precambrian  provinces  sug¬ 
gests  that  travel-tunes  observed  within  these  regions  might  be  transportable,  within  reasonable 
error,  to  other  similar  geological  provinces.  To  test  this  idea,  we  first  constructed  a  proposed  ^o- 
bal  Pn  reference  travel-time  curve  for  Precambrian  platforms  by  smoothing  and  averaging  die 
i^propiiate  Pn  travel-time  curves  of  Figure  34.  We  then  con^ared  this  proposed  global  Pn 
travel-time  curve  with  observed  data  firom  the  Australian  shield. 

Hgure  35  (bottom)  shows  this  composite  reference  Pn  curve  for  ancient  platforms,  as 
well  as  boundaries  indicating  two  standard  deviations  from  it  (s.d.  -  ±  0.72  sec).  The  concept  to 
be  tested  is  that  almost  all  Pn  travel-times  for  Precamluian  geological  provinces  should  lie  widiin 
this  interval.  This  conjecture  is  encouraged  by  current  theories  of  globtd  tectonics  and  continental 
drift  proposing  that  the  ancient  shield  and  platftum  areas  of  the  earth  have  a  common  age  and  ori¬ 
gin  as  a  single  unit  in  Pangea. 

Observational  data  came  fix>m  long-range  DSS  profiles,  and  from  recordings  by  seismo- 
logical  networks  of  large  chemical  and  nuclear  explosions,  carried  out  within  Australia.  Results 
of  these  studies,  including  record  sections  and  travel-time  measurements,  have  been  publisl^  in 
numerous  papers  (Bolt  er  n/.,  1958;  Doyle  and  Everingham,  1964;  Denham  era/.,  1972;  Hnlayson 
et  al.^  1974;  Muirhead  and  Cleary,  1977).  These  data  (Hgure  35,  top,  dots)  were  used  to  produce 
a  composite  Pn  travel-time  curve  for  the  Australian  Shield  based  on  explosion  seismolo^  obser¬ 
vations  (Figure  35  top,  broken  line).  The  standard  deviation  of  travel-time  measurements  from 
this  regional  Pn  curve  for  the  Australian  shield  is  ±0.79  sec.  A  second  Pn  travel-time  curve  for 
the  Australian  shield,  based  on  analysis  of  earthquakes  whose  locations  and  origin  time  were  well 
determined  [Bowman  and  Kennett,  1993],  is  alro  shown  on  Hgure  35  (top,  dash-dot  line).  The 
Pn  curve  inferred  from  analysis  of  earthquake  recordings  is  approximately  1-3  sec  earlier  than  the 
explosion  seismology  observations,  but  generally  supports  the  explosion  curve. 

When  the  Australian  regional  curve  is  superimposed  on  our  proposed  global  Pn  reference 
curve  for  Precambrian  shields  (Figure  35,  bottom),  it  may  be  seen  to  lie  within  the  ±2  s.d.  bound¬ 
aries  over  its  entire  range  (The  Pn  travel-time  curve  based  on  Australian  earthquake  recordings  is 
also  located  within  the  ±2  s.d.  interval  up  to  2000  km  distance).  As  a  result  of  this  test,  we  suggest 
that  our  composite  global  Pn  travel-time  curve  can  be  taken  as  the  initial  {q)proximation  for  other 
Precambrian  provinces  where  well-established  regional  curves  do  not  exist  We  note  that  accord¬ 
ing  to  numerous  geological  studies,  the  Precambrian  geological  provinces  constitute  over  half  of 
the  continents  and  include  a  number  of  regions  where  proliferation  may  be  a  concern. 

Conclusions 

For  the  Baltic  shield  and  neighboring  regions,  analysis  of  DSS  profile  observations  and 
recordings  of  mine  blasts  have  enabled  us  to  derive  a  new  Pn  travel-time  model,  a  new  method  of 
determining  regional  and  local  path-dependent  travel-time  corrections,  and  improved  phase  iden¬ 
tification  for  regional  events,  which  taken  together  promise  to  improve  the  accuracy  of  seismic 
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Figure  3S:  Regional  Pn  travel-time  curves  for  the  Australian  shield  inferred  from  explosion 
and  earthquake  observations  (fop)  and  proposed  globally  averaged  composite  Pn 
travel-time  curve  for  ancient  platforms  (bottom).  Dots  on  the  t(^  figure  are  obser¬ 
vations  from  explosions,  and  were  used  to  define  the  curve  shown  as  a  dashed  line. 
A  separate  curve  based  on  earthquake  data  is  also  shown  for  comparison.  The  Aus¬ 
tralian  explosion  curve  fits  within  the  ±2.  s.d.  boundaries  on  each  side  of  the  pro¬ 
posed  Pn  travel-time  curve  for  Precambrian  shields,  as  shown  on  the  bottom 
figure. 
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locations  in  Scandinavia.  Results  of  tests  completed  on  IMS  data  from  Kinina  mine  bla^  loca¬ 
tions  incorporating  this  new  information  into  interpretation  of  ARCESS  data  ^wed  iaqnove- 
ments  in  the  accuracy  of  their  locations,  although  sigiuficant  azimuthal  bias  remained,  and  there 
was  still  considerable  scatter  resulting  from  variations  in  the  semi-automatic  phase  time  picks. 
Results  can  probably  be  improved  if  we  use  several  arrays  for  event  locations,  better  define  the  Sg 
and  Sn  travel-time  curves  and  include  PmP  as  well  as  SmS  in  the  location  procedures. 

Analysis  of  explosion  seismology  observations  from  different  large  geological  provinces 
of  the  world  showed  that  large  regional  travel-time  variations  of  the  order  of  10  sec  fix'  Pn-waves 
(and  probably  15-18  sec  for  5/t-waves)  can  be  expected.  Pn  travel-time  curves  f<x  the  Rrecam- 
brian  platforms,  which  constitute  a  very  large  part  of  continents,  are  characterized  by  the  smallest 
Pn  travel-times.  Tectonically  active  regions  are  characterized  by  substantially  greater  travel- 
times.  \^th  these  large  and  region-dependent  travel-time  variations,  it  r^rpears  that  it  will  be  diffi¬ 
cult  to  produce  an  accurate  regional  epicentral  location  without  constructing  accurate  re^onal 
travel-time  curves  and  path-dependent  corrections  for  each  seismic  station. 

Travel-time  curves  for  Pn  in  Precambrian  shields  and  platforms  were  found  to  be  remark¬ 
ably  consistent  with  one  another,  permitting  us  to  define  a  composite  refererrce  Pn  curve  fen’ 
ancient  platforms,  with  standard  deviation  of  ±0.72  seconds.  To  test  the  transportability  of  this 
curve  to  other  platforms,  we  compiled  data  from  the  Australian  shield  and  found  that  it  agreed 
with  our  reference  curve  within  two  standard  deviations,  encouraging  us  to  believe  that  this  refer¬ 
ence  curve  may  be  applied  generally  to  similar  Precambrian  provinces. 

For  other  regions  of  the  continents,  including  tectonically  active  regions,  arrival  times  of 
Pn  are  generally  later;  variations  from  region  to  region,  however,  arc  very  large.  At  the  time  of 
writing  this  paper,  we  have  found  no  method  of  transporting  travel-time  curves  from  one  such 
province  to  another.  In  the  absence  of  transportability  principles,  it  appears  that  individual 
region-by-region  calibrations  will  be  required  to  obtain  good  location  accuracies. 
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css  Ground-TVuth  Database: 

Update  to  Version  1  Handbook 

L.  Grant,  J.  Coyne,  F.  Ryall 

Ground-Truth  Database  Overview 

This  report  is  an  update  to  the  CSS  Ground  Truth  Database  as  described  by  Grant  et  aL 
(1993).  The  objective  of  the  Ground-Truth  Database  (GTDB)  is  to  provide  events  of  known  type 
to  facilitate  research  on  the  seismic  discrimination  problem.  The  GTDB  is  a  collection  of  re^onal 
waveforms  and  carefully-reviewed  phase  parameters  generated  by  seimic  events  of  known  ^pe 
(Le.  earthquake,  quarry  blast,  etc.)  where  the  event  type  is  confirmed  by  some  means  in  additirm 
to  seismic  observations.  The  event  type  is  stored  in  the  etype  field  of  the  origin  table  and  is 
defined  as  follows. 


etype 

DESCRIPTION 

qb 

(^any  or  mine  blast  confirmed  by  quarry 

qb+ 

Quarry  or  mine  blast  with  designed  shot  inftnmation-  ripple  fired 

qb++ 

Quarry  or  mine  blast  with  observed  shot  information-  ripple  fired 

qbx 

(Quarry  or  mine  blast  -  single  shot 

qmt 

Quarry  or  mining-induced  events:  trenxrrs  and  rockbursts 

ex 

Explosion  of  known  origin;  e.g.  calibration,  cmistruction,  etc. 

nu 

Nuclear  explosion 

nc 

Nuclear  cavity  collapse 

eq 

Earthquake 

eq+ 

Earthquakes  in  a  swarm  or  aftershock  sequence 

eq++ 

Felt  earthquake 

o 

Other  source  of  known  origin 

u 

Undetermined  or  conflicting  information 

*  Note:  This  is  a  minor  revision  of  Tstbie  5  of  "CSS  Ground  Truth  Database”  \^:rsion  1  Handbook  (Grant  et 
al.,  1993)  (The  ex  definition  has  been  clarified.) 
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In  addition  to  the  above  changes  to  the  range  of  the  etype  field,  another  modification  to 
CSS  3.0  tables  is  that  the  ori^n  Table  is  used  to  sttxe  event  infotmatimi  collected  firom  experts  in 
the  local  areas.  This  has  resulted  in  a  "hybrid"  origin  Table  that  contains,  in  each  field,  die  best 
infoitnation  available.  To  indicate  this  mixed  origin  entry,  the  author  field  of  the  ori^  table  is 
set  to  "HYBRID". 

Currently,  the  GTDB  includes  288  known  events  located  in  Central  Europe  and  Scandina¬ 
via  and  recorded  on  stations  in  the  IMS2  and  GSETT-2  networks.  The  events  are  organized  into 
10  distinct  datasets,  based  on  geographical  location  and/or  source  of  informatitm. 

The  GTDB  is  available  through  the  ORACLE  database  account  called  "gtdb"  and  is  acces¬ 
sible  through  SQL^Plus  on  hugo  or  through  CenterView.  Version  1  and  Version  2  have  been  com¬ 
bined  into  the  gtdb  account  The  old  discriml  account  will  be  removed  on  July  1, 1994. 

Version  1  (Datasets  1-3),  previously  discriml  account; 

Version  2  (Datasets  4-10). 

The  GTDB  has  benefited  greatly  from  the  input  and  feedback  of  the  users.  This  project 
was  cancelled  in  January  of  1994  and  no  new  datasets  are  currently  being  added.  However,  please 
continue  to  notify  us  of  problems  encountered  in  using  the  data. 
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Version  1  Documentation 


The  documentation  for  Version  1  was  distributed  in  late  August  1993.  Ifyou  did  not  recieve  a  copy 
and  you  would  like  one  please  notify  Lori  Grant  (grant@seismo.css.gov  or  703-276-79(X)). 

title:  CSS  Ground  Truth  Database:  Version  1  Handbook 

number  Technical  Report  C}93-05 

date:  August  1993 

summary:  The  handbook  details  the  puipose  of  die  GTDB,  the  method 

of  building  it,  the  GTDB  Schema,  and  specifics  about  each  of 
the  three  datasets  comprising  Version  1.  It  also  contains 
a  plot  of  sample  waveforms  from  each  of  the  82  events. 


Version  1  Updates 

Two  of  the  datasets  in  Version  1  have  been  updated  since  publication  of  the  Version  1  Handbook. 
The  two  changes  briefly  noted  here  arealso  documented  in  the  log  file  in 

~grant/discriml_updates/discrim.log  on  the  machine  SOL. 

Dataset  #1 :  Vogtland 

The  etype  of  the  Vogtland  quarry  blasts  has  been  changed  from  qb 
to  qb+  to  indicate  that  there  is  some  additional  infcmnation 
available.  In  this  case,  the  total  charge  size  in  kg  is  stored  in  the 
minex  table.  The  minex  and  ori^n  tables  are  linked  through  evid. 

Dataset  #3:  Lubin 

New  locations  have  replaced  the  locations  published  in  the  Version 
1  Handbook.  All  location  calculations  have  been  poformed  by  the 
mines  on  the  basis  of  arrival  times  and/or  other  data  they  might  use 
(for  instance  the  location  of  damage,  if  any).  The  Institute  of 
Geophysics  converted  the  location  data  from  mine  coordinates  to 
geographical  coordinates.  The  error  in  previously  published 
locations  occurred  during  the  conversion.  The  Institute  of 
Geophysics  has  calculated  the  ori^n  times  and  magnitudes. 
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Version  2 


Overview 

Version  2  of  the  GTDB  comprises  7  datasets  fora  total  of 206 known  events.  10  events  of  unloown 
origin  are  also  included  because  of  diier  tenc^xxal  proximity  to  the  known  events.  die 
exception  of  die  Dataset  4,  all  resulted  by  obtaining  a  list  of  events  6om  a  published  paper  or  fitom 
the  author  of  a  published  paper  who  us^  the  events  in  a  study.  Widi  die  excqidon  ctf  Dataset  6 
the  events  within  a  given  dataset  are  actually  a  cluster  of  events  from  one  area.  Dataset  6  is  a  subset 
of  the  database  of  known  events  being  complied  at  the  Ruhr  University,  Bochum  in  ndiidi  events 
are  azimuthally  distributed  around  GERESS  and  generally  have  only  one  represteitative  event  in 
each  area.  Datasets  8-10  (121  events)  result  from  NORSAR  technical  reports  and  are  from  die 
^latity  mining  area  in  the  Khibinny  Massif. 


Waveform  data  are  from  the  IMS2  network  and  are  supplemented  by  some  GSETT-2  stations  for 
datsets  4  and  6.  Additionaly,  data  from  the  Apatity  3-C  station  was  supplied  by  NORSAR  staff  for 
the  quarry  blasts  in  dataret  #9. 

Starting  with  the  IMS2  bulletin,  parameter  data  for  each  event  has  been  carefully  re-ieviewed  by 
Flori  Ryall. 


Version  2  Datasets 

Dataset  #4 
number  of  events: 
event  type: 
evid  range: 
reference  papen 


Silesia _ 

31 

mining-induced  tremors  (qmt) 

373-1008 

none;  personal  communication  from  Institute 
of  Geophysics,  Polish  Academy  of  Sciences  in  die 
f(»m  of  a  list  of  all  events  >  ml  2.0  occuring  in 
Upper  Silesia  between  1  Jan  1991  and  30  June  1991. 


Dataset  #S 
number  of  events: 
event  type: 
cvid  range: 
reference  paper 


Swiss  Swanns . . 

27 

earthquake  swarms  in  Switzerland  and  Liechtenstein  (eq4-) 

1100-1142 

refid:  239 

author.  Kradolfer,  U. 

title:  RECENT  SEISMICITY  IN  SWITZERLAND  AS  LOCATED 
BY  THE  SWISS  NETWORK  AND  THE  GERESS  ARRAY 
yean  1992 

pub:  presented  at  GERESS  Symposium,  Waldkirchen,  Germany. 
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Dataset  #6 
description: 
numbtf  of  events: 
event  type: 

evid  range: 
reference  paper 


Dataset#? 
number  of  events: 
event  type: 
evid: 

reference  paper 


Gestennann . . . . . 

subset  of  events  from  Ruhr  University's  database 
14 

earthquakes,  quarry  blasts  and  mining-induced  tremors 
eq(4),  eq++(6),  qb(l),  qb+{l),  qiitt(2) 

1151-1169 
refid:  233 

author  Gestennann,  N.,  R-P.  Haijes,  M.  lost, 

J.  Schweitzer,  and  J.  Wuester 

title:  GERESS:  MONITORING  NATURAL  AND  ARTEFIOAL 
SEISMICITV  IN  CENTRAL  EUROPE ' 
year:  1992 

pub:  poster  presented  at  14th  PL/ARPA  Symposium,  Tuscon,  AZ 

Ruhr  Basin . 

13 

mining-induced  tremors  (qmt) 

1170-1184 
refid:  237 

author  Joswig,  ML  and  R  Schulte-Theis 
title:  MASTER-EVENT  CORRELTATONS  OF  WEAK  LOCAL 
EARTHQUAKES  BY  DYNAMIC  WAVEFORM  MATCHING 
year:  1993 

pub:  Geophys.  Joum.  Int,  vol  1 13 


Dataset  #8 
number  of  events: 
event  type: 
evid  range: 
reference  paper 


Apatity  Tremors . 

10 

mining-induced  tremors  (qmt) 

2000-2014 
refid:  240 

author:  Kremenetskaya,  E.  O.  and  V.  M.  Tijapitsin,  1992. 
tide:  INDUCED  SEISMICITY  IN  THE  KHIBINY 
MASSIF  (KOLA  PENINSULA) 
year  1992 

pub:  in  NORSAR  Scientific  Report  1-92/93 


Dataset  #9 
description: 

number  of  events: 
event  type: 
evid  range: 
reference  paper 


Apatity  Blasts-  3C . 

documented  quarry  blasts  recorded  after  apatity  3C  station 
upgraded  but  before  apauty  array  installed 
53 

quarry  blasts  with  charge  size  (qb+) 

2015-2067 
refid:  241 
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author  Mykkeltveiu  S. 

tide:  MINING  EXPLOSIONS  IN  raEKHmiNY  MASSIF 
(KOLA  PENINSULA  OF  RUSSIA)  RECX5RDED  AT  THE 
APATITY  THREE-COMPONENT  STATION 
pub;  Norsar  Scientific  Report  No.  14 


Dataset  #10 
description: 
number  of  events: 
event  type: 
evid  range: 
reference  paper 


Apatity  Blasts-  Array . 

documented  quarry  blasts  recorded  after  apatity  array  installed 
58 

quarry  blasts  with  charge  size  (qb+),  one  tremor,  (qmt) 

2070-2127 

refid:  244 

author  Kvaema,  T. 

tide:  INTELLIGENT  POST-PROCESSING  OF  SEIMIC 
EVENTS-PART2 

pub:  NORSAR  Scientific  Report  2-92/93 
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Explosion  Experiment  in  lymyauz 
August,  1993 

MYa.  Balbachan  and  M3.  Gokhbog 

E3q>eriiiient  Objectives 

Detecting  and  leccnding  the  mechanical-electric-magnetic  phenomena  (MEMP) 
that  occurs  within  rocks  as  a  result  of  destructitHi  by  an  explosion  are  the  main  objectives 
of  this  experiment  At  present  information  on  M^MP  caused  by  explositms  is  very  ccm- 
tradictoiy.  There  is  no  existing  reliable  data  on  the  petrc^hysical  ctxtditions  of  origin, 
source  and  physical  nature  of  MEMP.  It  was  necessary  to  carry  out  this  experiment  for 
MEMP  detection  at  different  distances  €r(»n  a  shot  point  in  the  attempt  to  obtain  such  data. 

Pressure  gradients  and  fractures  occur  in  the  rock  volume  simultaneously  with  an 
explosion  generating  seimiic  and  probably  electro-magnetic  effects  such  as  piezomag- 
netic,  electrcridnetic  and  so  on.  We  used  a  set  of  new  methods  developed  to  study  MEMP, 
to  detect  signals  and  to  get  infoimation  on  the  physical  nature  of  the  signals.  Electro-  and 
magnetotelluric  measurements,  near-surface  electric  measurements  and  seianic-acoustic 
observations  were  made  simultaneously  in  real  time  over  a  frequency  range  of  1  to  1000 
Hz. 


The  following  conditicms  were  fulfilled  for  massive  quarry  blasts  of  chemical 
explosions  having  yields  of  several  tons  to  several  tens  of  tons: 

1  •  source  dimensimis  and  energy  release  approximate  crustal  earthquakes  widi 

magnitudes  of  about  2.6  to  2.9; 

2  -  existence  of  shear  ccniqxrnent  essential  for  destructive  stress; 

3  -  simultaneous  generation  of  fractures  and  rock  crushing  within  the  focus  of  die 

explosions; 

4  -  splitting  of  electtic  charges  between  rocks  and  gas-dust  products  of  the  explo¬ 

sion. 

Electro-mechanical  emission,  electric  relaxation  due  to  fracture  generatirm  and 
friction,  seismoelectric  and  piezomagnetic  phenomena,  displacemmt  of  a  conducting 
body  in  the  geomagnetic  field,  variations  of  conductivity,  dielectric  penneability  and  other 
rock  properties  are  all  possible  sources  of  MEMP  discussed  in  scientific  literature.  Sepa¬ 
rating  electric  charges  between  the  rocks  and  the  gas-dust  cloud  generates  an  additional 
inmrfering  electrcxnagnetic  field.  This  kind  of  field,  generated  by  surface  explositms,  was 
identified  in  this  expoiment  by  using  specific  data  processing  and  was  not  observed  as  a 
result  of  underground  explositms. 
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Short  Review  of  Published  Data  on  Electro-Magnetic  Phenomena  Gener¬ 
ated  by  Explosions  in  Rocks 

This  is  a  short  review  of  the  electro-magnetic  phenomena  generated  by  explosions. 
Many  authors  have  observed  gas-dust  venting  generated  by  explosions,  accompMied  by 
electro-magnetic  fields.  According  to  V.N.  Shuleildn,  a  gas-dust  cloud  usually  has  a  posi¬ 
tive  electric  charge  and  the  rock  has  a  negative  one.  This  conclusion  has  been  confirmed 
by  remits  obtained  in  laboratories  by  other  authors.  Electric  potentials,  E,  of  iqrproxi- 
mately  10^  VAn  were  observed  within  vented  ckuids  generated  by  explosions.  Tl^  con¬ 
firms  the  possibility  of  recoding  E  in  air  at  distaiwes  varying  from  250  to  500  m  and 
measuring  origin  time  with  error  10*^  -lO*^  sec.  If  we  assume  that  the  ratio  of  dielectric 
permeability  in  air  (Ea)  and  in  soil  (Eg)  is  emial  to  the  conductivity  ratio  for  frequencies  of 
several  Hz,  in  other  words  is  equal  to  10*  ,  we  can  ignore  the  influence  of  the  electric 
field  generated  by  the  dust  cm  measurements  of  potential  differences  in  the  soil.  We  used 
this  coiclusion  fo  developing  the  experimental  technique. 

Electronagnetic  pulses  (EMP),  which  are  not  caused  by  discharges  in  the  air,  have 
been  recorded  by  some  authors  earlier  than  the  arrival  of  die  seismic  waves.  V.N.  Demin 
et  al.  (Sobolev,  Demin  et  al.,  1980)recorded  EMP  generated  by  underground  explosions  at 
polymetallic  me  deposits  within  the  frequency  range  of  0  to  10^  Hz  and  greater.  N.G. 
Khatiashvili  and  other  researchers  report  on  EMP  observations  between  10^  -10^  Hz 
caused  by  the  destruction  of  non-metadliferous  rocks  (Khatiashvili,  1989). 

M.B.  Gokhberg  et  aL  published  the  results  of  MEMP  observatimis  during  and  after 
seismic  wave  propagation  (Gokhberg  et  al.,  1987).  Smoetimes  it  can  be  explained  by  sen¬ 
sor  vibrations  relative  to  stationary  natural  fields.  Electric  waves  caused  by  seismic-elec¬ 
tric  pheiKxnena  will  be  genoated  simultaneously  widi  sdsmic  waves.  According  to 
Ya.M.  Frenkel  et  al.,  a  linear  relationship  should  be  observed  betweoi  amplitudes  of  the 
electric  and  seismic  waves.  Electro-magnetic  waves  generated  by  Raylei^  waves  (Gor¬ 
bachev  et  al.,  1987)  or  having  iKm-hartrxHiic  waveforms  and  low  frequency  1-10  Hz  have 
also  been  observed.  Results  of  observations  of  anomalies  of  quasistatic  magnetic  fields 
observed  during  1-10  days  after  explosimis  woe  published  (Oganesyan,  1986;  Belekon, 
1992;  Balbachan,  1988).  These  anonudies  can  be  generated  by  a  re^ual  piezomagnetic 
effect  Local  variations  of  electrocmiductivity  or  collector  properties  of  rocks  can  cause 
anomalies  of  the  electric  field  after  explosicms  (Adushldn  et  al.,  1989). 

This  review  shows  that  the  results  of  the  observations  of  electromagnetic  effects 
of  explosimis  are  contradictory  and  unreliable.  New  experiments  aimed  at  collecting  data 
mi  eiqilosimi  parameters,  on  electromagnetic  signals  generated  by  the  explosimis  and 
recorcted  at  different  distances,  and  on  characteristics  of  rocks  are  essential  to  resolve 
these  discrepancies. 


78 


The  Experiment  Conditions 

The  following  sciaidsts  participated  in  the  experinoent:  M.Ya.  Balbachan,  senior 
scientist,  was  in  charge  of  the  observation  technique,  nuuiagement  and  operation  oi  seis¬ 
mic  and  electrcwoetric  sensors;  S.M.  Krylov,  leading  scientist,  was  in  cha^  die  opera¬ 
tion  of  magnetcHneters;  and  I£.  Matveev,  senior  scientist,  was  in  charge  ci  die  recording 
equipment 

The  experiment  todc  place  frmn  19  to  31  August  1993  at  a  high-altitude  c^ien  pit  oi 
the  Tymyauz  tungsten-molybdenum  integrated  plant  near  the  lymauz  sttdement  Mtuiag- 
ers  and  enqiloyees  of  the  integrated  plant  khuUy  helped  us  to  fulfil  die  eiqperiment 
Between  22  and  29  August  two  explosi(»s  were  detonated.  One  of  them  was  an  under¬ 
ground  explosion  and  the  second  was  a  surface  explosicMt  The  underground  explosion 
was  detmiated  in  a  mine  immediately  beneath  the  quarry  at  9:00  a.m.  Moscow  time  and 
the  surface  explosion  was  detonated  in  the  upper  part  of  the  quarry  in  overburden  rocks  at 
10:00  a.m.  Moscow  time.  Receding  statimis  were  located  near  a  mountain  tc^  at  heists 
of  26SS-266S  m  above  sea  leveL  At  the  observation  area  there  was  fog  following  rain  and 
a  light  southerly  wind  during  the  first  explositm.  The  weather  was  dry  and  fine  during  the 
second  explosion. 

This  experiment  was  completed  within  a  rock  massif  c(»sisting  of  marble,  lime¬ 
stone  and  gamet-pyroxenite  with  molybdenite,  and  biotite,  pyroxene,  and  plagioclase  with 
quartz  hmnfels.  Small  druses  of  galena  crystals  were  observed  c«  marble  surfaces  within 
the  area  cf  measurements.  Rough  estimates  of  the  electric  conductivity  in  the  near-suri&ce 
layer  within  the  observation  area,  based  cm  short-base  measurements,  yielded  10*^  Sm  and 
10*^  Sm  for  humid  and  dry  weather,  respectively.  The  content  of  ferromagnetics  in  the 
rock  is  very  small,  and  iron  oxides  are  present  in  the  form  of  limonite  (Chief  geologist 
pers.  ccxtl).  The  rock  massif  should  be  paramagnetic  in  general.  Rock  density  is  2.87  g/ 
cm^  on  average. 

Both  explosions  were  detonated  on  Sunday  when  normal  work  at  the  underground 
and  surface  pits  was  not  performed.  All  power  sources  in  the  mine  and  pit  were  switched 
off  during  tire  explosion  recordings.  Ihe  closest  point  where  electric  power  was  not 
switched  off  was  the  pit  office  building  located  at  a  (fistance  of  approximately  one  kra  A 
cable  ore-conveyor  was  discotmected  during  the  explosions.  The  top  station  of  the  cable 
conveycm  is  legated  near  the  pit  office  building.  There  were  no  unambiguous  signs  of 
techrucal  noise  in  the  specific  reccredings  made  to  study  noise. 

An  underground  explosion  was  ctetonated  on  August  22  at  three  points  simulta¬ 
neously  within  levels  between  2150  m  and  2410  m  above  sea  level.  The  stnallest  distance 
of  the  shot  point  fiom  a  recording  station  was  approximately  450  m.  The  g«ieral  yield  of 
the  explosion  was  greater  than  10  tons.  Because  of  techrucal  reasons,  this  explosion  was 
notreoitded. 

A  surface  explosion  was  also  detonated  on  August  22  at  one  point  fiv  stripping. 
Explosive  charges  of  350  kg  of  grammonite  were  detonated  in  drilled  holes  in  two  heni- 


zcms  at  2665  m  and  2655  m  above  sea  level,  and  at  depths  of  tqpproximately  10  m.  Drilled 
holes  were  located  in  11  rows.  The  distance  between  the  rows  and  the  drill  hrdes  within 
the  rows  was  3  rtL  Drilled  holes  in  each  row  were  detonated  simultaneously,  and  betwera 
the  rows  there  was  a  20  msec  delay.  The  total  explosion  time  was  240  msec  and  the  total 
yield  was  approximately  32  ttms.  The  total  vtdume  of  rocks  extracted  was  13000-15000 
m^  and  the  ^struction  zone  volume  was  35x15x10  m  -i-  39x10X10  m.  The  total  area  of 
fractures  generated  by  the  explosion  is  estimated  at  ^proximately  10,000  -  100,000  m^. 
The  distance  from  the  shot  point  to  the  recording  station  was  350  m. 

An  urulerground  explosion  was  detonated  at  a  depth  of  100  m  from  the  surface  on 
August  29  at  a  point  located  at  2302  m  above  sea  level.  The  distance  between  the  shot 
point  and  the  recxirding  stations  varied  from  600  m  to  650  m.  The  yield  of  the  explosion 
was  5.8  ums.  The  total  volume  of  rocks  extracted  was  3,000  m^,  and  the  total  area  of  frac¬ 
tures  generated  by  the  explosion  was  approximately  2,000-20,000  m^.  The  recording  sta¬ 
tions  woe  locat^  along  a  line  at  2665  m.  This  line  was  approximately  normal  to  a  line 
cotmecting  the  shot  point  and  the  first  (northern)  recording  station. 

The  surface  explosion  detonated  on  August  29  was  located  at  2675  m.  The  charges 
were  distributed  simil^y  to  the  explosion  demnated  on  August  22  and  the  yield  was  7.3 
tons.  The  total  explosion  time  was  100  msec.  Destruction  zone  volume  was  approxi¬ 
mately  20x15x10  m.  Total  volume  of  rocks  extracted  was  3,000  m^.  The  total  area  of 
fractures  generated  by  the  explosion  was  approximately  2,000-20,000  vc?.  The  distance 
between  the  shot  point  and  the  first  recording  station  was  250  m.  The  shot  points  were 
located  at  the  same  points  as  the  underground  explosion  and  at  the  same  radius  to  the 
south  from  the  shot  point  at  2665  m. 

The  objectives  of  these  field  observations  were  to  check  the  sensitivity,  the  current 
level  of  the  recording  system,  and  the  existence  of  electrical  and  magnetic  effects  gener¬ 
ated  by  destruction  of  the  rock  by  the  explosions.The  electric  and  magnetic  sigruds  can  be 
recorded  earlier,  simultaneously  and  later  than  seismic  signals.  The  study  of  die  space 
attenuation  of  the  anticipated  electric  and  magnetic  signals  was  another  aim  of  this  experi¬ 
ment 

Methods  of  Measurements  and  Parameters  of  Recording  Equipment 

The  recording  stations  were  located  along  the  radius  from  the  center  of  the  surface 
explosions.  The  distance  between  the  points  of  tire  recordings  for  die  explosion  on  August 
29  was  1(X)  m  and  the  distances  between  the  shot  point  and  the  recording  stations  varied 
from  25()m  to  550  m  (Hgure  36  top).The  only  recording  station  for  the  explosion  on 
August  22  was  located  at  a  distance  of  350  m  from  the  center  of  the  surface  explosion. 
One-component-electric  and  one-component-magnetic  fields  were  measured  in  the  air  or 
in  the  soil. 

Each  recording  station  included  a  magnetometer  designed  by  S.M.  Krylov.  The 
magnetometer  was  installed  on  August  22  horizontally  and  oriented  towards  the  magnetic 
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pole.  Three  similar  magnetometers  were  also  installed  for  recording  the  August  29  explo¬ 
sion  at  points  1-3  and  were  oriented  along  a  major  geomagnetic  vector.  In  both  cases,  hor¬ 
izontal  projectimis  of  the  magnetometer  axis  were  oriented  tqiproximately  alrnig  the  radius 
towards  die  surface  explosicm  centers  and  popendicularly  to  the  center  of  the  underground 
explosions.  Spherical  ctyiacitance  antennas  were  installed  3m  from  the  magnettHoeters  at 
1  m  height  above  the  Earth’s  surface  to  measure  the  vertical  ctxnponent  of  Ea  and  to  mea¬ 
sure  the  radial  component  of  the  horizcmtal  Ea  with  help  of  two  similar  antennas. 

Accelerometers  were  installed  immediately  beneath  the  antennas  on  solid  hard 
rock.  A  60°  angle  was  produced  in  the  vertical  plane  between  the  measuring  axis  of  the 
magnetometers  and  a  radus  to  the  surface  explosion.  A  measuring  electrode  was  buried  at 
a  distance  of  1  m  from  each  accelerometer.  These  electrodes  formed  dipoles  to  measure 
the  difference  of  potentials  Ug  at  points  1, 2, 3  and  4  (Figure  36  bottom).  Spherical  anten¬ 
nas  and  measuring  electrodes  were  cormected  with  similar  amplifiers  having  3x10*^  Ohm 
input  resistance.  Signals  from  magnetometers,  accelerometers  and  amplifiers  at  points  1 
and  2  fed  input  of  the  blocks  consisting  of  four  charmels.  The  blocks  performed  the  fre¬ 
quency  modulation  of  signals  and  generated  a  pilot  signal  to  make  a  time  scale  and  to 
eliminate  errors  caused  by  uncontrolled  variations  of  recording  velocity.  Ouq)ut  signals 
firom  the  magnetometer,  the  accelerometers  and  the  amplifiers  fed  input  to  a  TEAC  XR- 
3()C  8-chatmel  magnetometer  for  point  3  and  point  4.  Sound  accompaniment  was 
recorded  by  charmel  8. 

An  ultrashort-  wave  radio-synchronizer  and  radio-receivers  that  were  cormected  to 
a  tape-recOTder  were  used  to  syncluronize  the  recordings  from  different  points.  The  syn¬ 
chronization  pulse  was  recorded  several  seconds  earlier  dian  the  surface  explosion  origin 
time.  This  pulse  and  the  "pilot  signal"  provided  rather  accurate  information  on  origin 
time.  We  assumed  that  we  would  be  able  to  reccxd  a  signal  generated  by  the  gas-dust 
venting  with  the  help  of  the  spherical  antetma.  Two  12-volt  alkaline  batteries  were  used  as 
the  power  source  for  each  point  of  the  recording.  Grounding  was  made  for  each  point  3  m 
firom  the  measuring  electr^es  using  a  similar  electrode  and  increasing  the  NaQ  concen¬ 
tration.  All  of  the  equipment  was  isolated  firom  the  Earth’s  surface. 

Field  observations:  All  elements  of  the  recording  system  woe  installed  sevoal 
hours  before  the  explosion.  The  noise  was  recorded  during  the  O.S-1.0  hour  after  the 
power  and  all  of  the  equipment  at  the  pit  and  the  mine  were  cut  off.  The  tape-recorders 
were  turned  on  5  min.  before  the  tinre  of  the  explosions  and  were  turned  off  15  min.  later. 
Recording  equipment  was  calibrated  afto*  the  explosion  recordings  by  signals  generated 
as  a  result  of  10  kg  weight  drops. 

Equipment  parameters; 

1.  Magnetometer.  Signal-to-noise  ratio  for  1  Hz  was  greater  than  7x10^.  The  peak 
amplitude  value  of  the  output  signal  was  ±8  V. 

2.  Spherical  antenna.  Capacitance  with  an  unscreened  connecting  wire  was  IS  pF. 
The  height  above  the  Earth’s  surface  was  1  m. 
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3.  Accderometer.  The  nmninal  conversion  coefficient  was  0.S  V  sec^/mu  The 
maximal  signalAioise  ratio  was  not  less  than  106  dB.  The  peak  amplitude  value  of  the  out¬ 
put  signal  was  ±10  V.  The  ouq)ut  resistance  was  2  kOhm. 

4.  The  measuring  electrodes  were  mack  firran  indium  and  had  the  shape  of  a  cross. 
The  electrode  parameters  woe  100x100x100  mm. 

5.  Block  consisting  of  two  aiq>lifiers.  Input  resistance  was  3  MObm.  The  coeffi¬ 
cient  of  suppressing  for  SO  Hz  agnals  was  10^.  These  anqplifiers  were  idmtical  and  inde¬ 
pendent  The  maximal  coefficient  of  the  amplificatitm  was  65  Db.  The  maximal  ouq>ut 
signal  was  ±S  V.  The  maximal  signal-to-noise  nuio  was  not  less  than  SxlO^. 

6.  The  differential  preliminaiy  amplifiers  were  included  between  the  amplifiers 
mentioned  above  (see  point  S.  Figure  36  bottcnn)  and  the  electrodes.  Input  resistance  was 
860  kOhm.  The  maximal  ou^ut  signal  was  ±10  V.  The  maximal  signal-to-noise  ratio  was 
not  less  than  10^. 

7.  Four-channel  blocks.  The  dynamic  range  for  the  firequency  noodulatitm  was  40 
dB.  The  central  frequency  of  the  modulator  was  7.5  Idiz.  The  pilot  signal  frequency  was 
2.5  kHz.  The  maximal  ouqtut  signal  was  5V. 

8.  Magnetograph.  The  parallel  recording  of  8  channels  to  video-cassette  was 
completed.  The  recording  rates  were  2.4  sm/sec  and  4.8  sm/sec,  respectively.  Frequency 
modulatitm  was  used.  The  dynamic  range  was  45  dB.  Frequency  range  varied  firm  0  to 
625  kHz  and  from  0  to  1250  kHz  for  recording  rates  of  2.4  Wsec  and  4.8  sm/sec,  respec¬ 
tively.  Maxiroal  ouqtut  signal  was  1 V. 

9.  Tfqte-recordets  were  used  with  two  channel  parallel  recording  to  an  audio-cas¬ 
sette.  Dynamic  range  was  40  dB  fOT  frequency  7.5  kHz.  Maximal  input  signal  was  5  V. 

10.  Radio-synchronizer  included  a  harmonic  signal  transmitter  for  a  frequency  of 
KXlMHz.  Its  power  was  2.5  Wt  The  synchronic  pulse  was  generated  by  a  sharp  cut  off. 

These  parameters  of  the  recording  equipment  archived  the  following  sensitivity  of 
measurements:  acceleration  of  seismic  signals  -  KT^  m/sec^;  magnetic  field  (magneto- 
gnq>h  recordings)  -  Bi=2KT  for  frequency  1  Hz  and  Bis2nT  for  frequency  16  Hz;  die  elec¬ 
tric  field  in  air  for  frequency  1  Hz  Eai^.l3  Vfrn  for  magnetograph  recordings  and  1.3  V/ 
m  for  tape-recorder  recording;  electric  field  in  air  for  frequency  KXX)  Hz:  1.3x10*^  VAn 
for  magnetogr^h  recordiira  and  1.3xl(r^  V/m  for  t^-recoider  recording;  for  an  electric 
field  in  a  soil  Eods2.5xlO'^V/m. 

Noise  Measurements 

Noise  was  recorded  continuously  by  the  magnetognqihs  during  the  two  hours 
before,  after  and  between  the  explosions  detonated  on  August  29,  1993.  The  magneto- 
grqih  has  <»ie  order  greater  sensitivity  in  comparisrai  with  a  tape-reconder  and  we  used  the 
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magnetograph  recwdings  of  noise  for  analysis.  As  we  mentioned  above,  sources  of  ncnse 
generated  by  human  activity  were  eliminated  within  a  distance  greater  than  1  km  fipom 
recording  points.  Dipoles  at  point  3  and  pc^t  4  had  a  100  m  base  and  were  aligned  fiom 
north  to  south.  The  similar  con^nent  oi  the  electric  field  in  the  soil,  the  signals  fiom  two 
^helical  antennas,  the  two  accelenxneters  and  one  magnetometer  were  recorded  (point  3, 
Figure  36  bottom). 


'IMth  acceleixHiieters  we  were  able  to  record  the  seismic  signals  generated  by 
explosions,  by  the  pec^le  near  the  recording  points  and  by  the  fall  of  the  sttmes  after  tiw 
exploskms.  ^th  the  exception  of  these  nxmients,  noise  signals  were  manifested  by  die 
internal  noise  of  the  equipment  Signals  recorded  by  other  sensors  were  very  ccmiplicated. 

1.  High-frequency  continuous  signal.  This  signal  was  recorded  bv  all  of  the  elec- 
trcHnagnetic  sensors.  The  dcxiiinant  frequency  of  signals  recorded  with  the  help  of  the 
electrodes  and  antennas  was  approximately  300  Hz.  Requency  qiectra  varied  slowly  with 
time,  important  variations  happened  during  minutes  or  tens  of  minutes.  The  typical  range 
of  the  amplitude  variations  was  observed  fimn  10*^  to  10*^  V/m  and  finom  2.5xl(r^  to  10'^ 
V/m  for  signals  recorded  with  the  help  of  the  spherical  antennas  and  electrodes,  corre¬ 
spondingly.  The  general  source  of  the  signals  located  at  distances  larger  than  100  m  man¬ 
ifests  by  synchronism  of  signals  receded  at  different  points  and  amplitude  mnilarity.  The 
signals  recorded  by  electrodes  have  a  90^  phase  s^  in  ccmqrarison  with  the  signals 
reaxded  with  the  help  of  ^herical  antennas  that  fit  a  phase  shift  between  the  current  of 
displacement  and  the  current  of  conductivity. 


2.  Harmonic  low-ftequencv  signal.  This  signal  was  recorded  by  all  of  the  electro¬ 
magnetic  sensors.  The  signals  recorded  by  the  magnetometer  and  electrodes  immediately 
before  the  underground  explosion  detoruUed  on  August  29  had  the  dominant  fiequency  of 
0.2  Hz  and  amplitudes  of  SO  xT  and  5x10'^  VAn,  ctxresptxidingly.  The  signals  recorded 
by  spherical  anteimas  had  fiequencies  of  0.067  Hz  and  1.0  Hz  and  anqrlitudes  of  0.5  VAn 
and  0.1  VAn,  correspcmdingly.  The  signals  recorded  by  electrodes  and  magnetometers 
during  44  minutes  before  the  surface  explosion,  had  a  finniuency  of  0.2  Hz  and  an^litudes 
of  100  Ttr  and  SxlO*^  VAn;  meanwhile,  signals  fiom  the  spherical  antennas  had  fiequen¬ 
cies  of  0.2  Hz  and  0.7  Hz  with  corresponding  anqrlitudes  of  0.5  V/m  and  0.2  VAn.  Hiase 
diift  between  signals  from  the  electrodes  and  the  antetmas  was  approximately  60^;  that 
can  be  explained  by  the  influence  of  a  conductivity  current  in  the  air  fen:  the  mentioned  fie¬ 
quencies. 


3.  "Steps"  for  the  magnetometer  recordings.  The  steps  had  a  positive  and  a  neea- 
tive  polarity.  Amplitude  reduction  with  time  was  exponential  with  the  relaxatiem  period 
equal  to  approximately  7  sec.  This  is  a  typical  value  of  internal  telaxafitm  poiod  for  the 
magnetometer.  The  amplitudes  of  the  steps  varied  fiom  10  to  50  xT.  The  typical  varia¬ 
tions  of  the  records  fiom  the  electrodes  and  the  antennas  corre^nded  to  the  steps.  Atnoo- 
spheric  signals  that  were  greater  than  the  noise  level  were  observed  on  the  average  tmee 
per  minute  for  the  spherical  anterma  and  the  magnetometer  rectndings.  Mtxiochromatic 
harmonic  80  Hz  signals  were  recorded  by  the  magnetometers  and  were  not  recorded  by 
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the  aocekrooieten. 


Eiplosion  Effects 

Electnmiagnetic  sigiuils  before  the  seismic  waves  were  not  reooided  for  any  of  die 
dnee  explosioiis  widun  die  dynamic  range  of  obsexvadoiis.  Anomalres  were  observed  on 
spherical  antenna  recordings  simultaneously  widi  die  venting  of  die  gas-dust  doud  for 
surfoce  explosions  ^gure  37).  Arrows  on  diis  and  other  figures  mark  the  origin  time  oi 
ex]riosions  and  die  arrival  time  of  seismic  waves  at  recording  points. 

The  following  effects  occurred  during  and  after  die  propagation  of  sdsmic  waves. 

1).  The  seismic  signal  caused  the  surface  eaqilosion  detonated  on  August  29, 
1993  had  a  waveform  typical  of  nmldple  evoits  uddi  die  second  eaqilosion  ddayed  by  1.4 
sec.  The  iduue  velocity  calculated  fitom  the  arrival  time  of  the  seismic  signal  measured  at 
pmnt  3  and  pmnt  4  was  4850  nVsec  (E^gnre  38  and  Hgure  39).  The  sdsmic  signal  gener¬ 
ated  by  the  underground  explosion  detonated  on  August  29,  1993  is  shown  in  Hgure  40 
and  41. 


2) .  The  magnetic  field  increased  sharply  during  seismic  wave  propagation  (Hg- 
ures  42-45  and  Hgures  46-47).  The  relaxation  time  observed  for  diere  anomalies  was 
tpprmdrruttdy  7  sec.  This  time  was  influenced  by  the  recording  system  and  we  can 
assume  that  teal  time  was  greater.  The  signal  recorded  by  die  magnetometer  from  die 
oqdosion  detonated  on  August  22  had  an  essentially  different  waveform  (Hgure  48)  in 
coopatison  with  die  explosion  detonated  on  August  29  (Hgure  42).  That  can  be 
explained  by  different  sensor  orientation.  The  attenuation  oi  die  sdsmic  wave  anplitude 
was  observed  vs.  distance.  The  results  of  calibration  showed  that  the  modulation  oi  the 
agnal  recorded  at  point  3  fitom  the  explodcm  detonated  of  August  29  was  caused  by 
mechanical  resonance  of  the  sensor. 

3) .  The  results  of  calibration  hqpened  to  be  very  interesting.  The  noagnetonaerers 
manifested  a  stqr-shipe  increase  of  die  noa^ietic  field  at  different  observation  points  simi¬ 
lar  to  the  eiploaondfecL  Avery  pecific  effect,  "memory"  was  observed  during  die  cali¬ 
bration  vdien  die  first  blows  only  generated  an  observed  jdienomena  characterizing  a  12 
sec  relaxatkm  period.  The  jdnt  analysis  ci  ndse  recordings  and  the  results  d  the  record¬ 
ing  system  tests  yielded  die  informatimi  dut  die  jdienomena  can  not  be  eiplained  by  an 
influence  of  die  recording  equipment  (Hgures  49  and  50) 

4) .  Amplitude  decreasing  of  background  signals  were  observed  simultaneously 
widi  the  previously  described  variations  of  the  nugnetic  field  measured  by  dectrodes  and 
qdierical  antennas.  This  phenomenon  was  observed  up  to  die  relaxation  of  die  magnetic 
field  signal  (Hgures  51-54  and  Hgures  55-58).  It  is  wordi  noting  diat  stqi-by-siep  pulses, 
B,  having  an  anplitude  equal  to  tens  of  ifUs  remained  during  the  anonaalies  generated  by 
die  explosions.  The  origin  time  of  the  explosions  detcxuued  on  August  28  coinckie  with 
die  time  arrivals  of  dieae  pulses  (Hgures  46  and  47). 
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5).  The  decrease  in  the  ratio  Eg/Ea  for  the  background  signal  of  low  fiequency 
during  the  44  nunutes  b^ween  explosiois  can  be  e3q>lained  by  die  explosion  influence. 

Discussion  of  Results 

Results  of  these  experiments  showed  that  widiin  measurement  errors,  and  for  die 
frequency  range  1-1000  electromagnetic  signals  generated  by  explosicms  and  arriving 
earlier  than  seismic  signals  were  not  observed.  This  result  corresponds  to  the  source 
radius  r^  ~  IS  m,  distance  from  die  explosion  center  Ro  =  3S0-6S0  m,  the  eiqiloritm  mag¬ 
nitude,  and  rock  conqiosition.  If  we  assume  that  the  source  is  an  electromagnetic  dipde, 
dun  is  the  simplest  pattern  widi  the  smalle;*t  attenuation  of  die  field  vs.  distance,  die  elec¬ 
tric  field  in  the  air  above  the  source  will  be  Eoo  »  0.S  Ecu  x.  (Rofro)^  <  1800  VAn  for  1 
Hz  and  18  mV/m  for  1  kHz,  where  Eoti  is  the  maximal  sensitivity  of  measuremrats  for 
underground  explosions.  The  typical  value  of  the  electric  field  at  die  source,  taking  into 
account  the  size  of  the  measuring  dipole,  is  Ego  =  Egi  (Rofro)^  <  9x10*^  V/m. 

Meanwhile,  Bo  =  Bi  (Ro/ro)^  <  1.8  xT.  Actually,  local  fields  can  have  far  larger 
values  but  attenuatitm  fields  vs.  distance  will  be  far  larger.  The  absence  of  electromag¬ 
netic  fields  recorded  earlier  than  seismic  wave  arrivals  shows  that  the  eiqilositm  source  is 
not  the  source  of  MEMP  with  far  greater  qiecific  power  than  rock  samples  destroyed  in 
laboratories  and  generating  usually  smaller  electmnagnetic  fields.  Electromagnetic  radia- 
ticm  can  be  receded  near  the  source  <xdy. 

We  believe  diat  events  observed  during  this  experiment  and  recorded  by  electro¬ 
magnetic  sensors  can  be  explained  in  terms  of  the  tenporary  increase  of  mountain  mass 
conductivity  caused  by  seismic  wave  jaxpagation.  As  restdt,  the  density  of  the  telluric 
current  and  the  seccxidary  magnetic  field  will  increase.  On  the  other  hand.  Eg  will  be 
decreased  and  high-fi«quency  signals  will  be  screened.  Anplitude  modulation  will  follow 
this  process.  The  sharp  increase  of  mountain  mass  ctmductivity  after  die  underground 
explosion  detonated  (m  August  29  is  not  strange. 

The  increasing  of  rock  conductivity,  probably  caused  by  mechanical  and  cheocdcal 
factors,  was  observed  after  vibrator  performance  Ofolkova  et  al.,  1988)  and  in  laboratories 
when  acoustic  waves  prcpagated  tb^gh  marble  sanples  (Sadovsldi  et  al.,  1989).  This 
phenmaencm  has  not  yet  been  explained  frcmi  the  point  view  of  physics.  Further  investi¬ 
gations  in  this  field  will  give  us  an  additional  chance  for  the  remote  study  of  seismic 
waves  within  the  near  zone  from  the  stnirce  and  the  understanding  of  MEMP  antnnalies 
observed  from  earthquakes. 
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Figure  36:  The  experiment  scheme  performed  on  August  29, 1993. 
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Figure  37:  Surface  explosion  detonated  on  August  22, 1993,  point  No.  1  located  at 
a  distance  of  3S0  m  from  the  explosion  center.  Variations  of  the  vertical 
component  of  the  electric  field  in  air  Ea  vs.  dme. 
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Hgure  38:  Sinface  explosi<m  detonated  on  August  29, 1993;  point  No.  3.  Variaticms 
of  seismic  wave  acceleruicm  vs.  time. 
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Hgine  39:  Surface  explosion  detonated  on  August  29, 1993;  point  No.  3.  Variations 
of  seismic  wave  acceleratimi  vs.  time.  High  frequency  filter  was  used. 


Figure  40:  Undeiground  explosi<m  det(Miated  on  August  29,  1993;  point  No.  3. 
Variations  of  seismic  wave  acceleration  vs.  time. 

DC4 

JQCQ 

0 

-102 

-JIM 


« 

0  1  2  3  4  5  6  7  5  9  10  11  12  13  14  15  K  17  18  19  30  21  22  23  24  • 

t,  s 


Ingure41;  Undeiground  explosion  detonated  <xi  August  29,  1993;  point  No.  3. 

Vuiations  of  seismic  wave  acceleraticm  vs.  time.  High  frequency  filter 
was  used. 
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Figure  43:  Surface  explosicm  detonated  on  August  29, 1993;  point  No.  2  located  at 
350  m  distance  from  the  explosion  center.  Variations  of  die  magnetic 
induction  B  vs.  time. 
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Hgure  44:  Surface  explosicm  detonated  on  August  29, 1993;  point  No.  3  located  at 
450  m  distance  from  the  explt^tm  center.  VariadtHis  of  the  magnetic 
induction  B  vs.  time. 
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Figure  46:  Undoground  explosion  dettmated  on  August  29,  1993;  pcwt  No.  3 
located  at  650  m  frcxn  the  explosion  center.  Variations  of  the  magnetic 
induction  B  vs.  time. 
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Figure  47:  Undeiground  explosion  detonated  cm  August  29,  1993;  point  No.  3 
kmatedat  650  m  from  the  explosion  center.  Variations  of  the  magnetic 
inducticm  B  vs.  time.  Low-frequency  filter  was  used. 
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Figure  48:  Surface  explosicm  detonated  chi  August  22, 1S>93,  point  No.  1 .  Vuiadons  of 
the  radial  component  of  the  magnetic  induction  B  vs.  time. 
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Hgure  49:  Calibrating  agnals  leomted  on  August  29,  1993  at  point  No.  2.  Varia 
tions  of  magnetic  inductkm. 
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Hgure  50:  Calibrating  signals  receded  on  August  29,  1993  at  point  No.  2.  Varia 
ti<xis  of  seismic  wave  acceleration. 
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Hguie  SI:  Surface  explosi<»  detcmated  on  August  29, 1993;  ^oint  No.  3.  Variations 
of  the  radial  om^onent  of  the  electric  field  in  soil  Eg  vs.  time. 
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Hgure  52:  Surface  explosion  detonated  on  August  29, 1993;  point  No.  3.  Variations 
of  the  radial  component  of  the  electric  field  in  soil  Eg  vs.  time.  Low-fre¬ 
quency  filter  was  used. 
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Figure  53:  Surface  explosion  detonated  on  August  29,  1993;  point  No.  3. 

Variations  of  the  vertical  component  of  the  electric  field  in  soil  Ea 
vs.  time. 
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Figure  54:  Surface  explosion  detonated  on  August  29,  1993;  point  No.  3. 

Variations  of  the  vertical  component  of  the  electric  field  in  soil  Ea 
vs.  time.  Low-frequency  filter  was  used. 
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Figure  55:  Undapound  explosioa  detcuuied  on  August  29, 1993;  p<wt  No.  3.  Van- 
ati<ms  of  the  radial  conqponent  of  die  electric  field  in  soil  Eg  vs.  time. 
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Figure  56:  Underground  explosion  detonated  on  August  29, 1993;  point  No.  3.  "Wm- 
atkms  ci  the  radial  conqx>nent  of  die  dectric  field  in  soil  Eg  vs.  time. 
Low-fiequency  filter  was  used. 


99 


Eo ,  V/m 


Figure  57:  Underground  explosion  detonated  on  August  29, 1993;  point  No.  3.  Vari¬ 
ations  of  the  vertical  con^nent  of  the  electric  field  in  air  vs.  time  Ea. 


Figure  58:  Underground  explosion  detonated  on  August  29, 1993;  point  No.  3.  Vari¬ 
ations  of  the  vertical  component  of  the  electric  field  in  air  vs.  time  Ea. 
Low-frequency  filter  was  used. 
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Filming  of  the  Ripple-Fired  Explosion  at  the 
lyrnyauz  Industrial  Complex 

By  V.L  'nykanovsidi 
Institute  for  Dynamics  of  the  Geo^heies 
Russian  Acadraiy  of  Sciences 


Abstract 

High  speed  filming  of  the  ripple-fired  explosion  at  the  open  quaixy  of  die  lymyanz 
industrial  complex  (North  (Caucasus)  has  been  carried  out  Objectives  of  tiiis  project  were 
to  obtain  information  on  the  detonation  of  a  sequence  of  charges  in  drilled  holes,  on  real¬ 
time  delays  and  other  specific  features  of  a  set  of  chemical  explosions. 

Project  description 

According  to  the  explosion  plan,  two  sets  of  grammonite  chemical  charges  of  total 
weight  4S  tons  were  to  be  detonated  at  the  lymyauz  wolfiramite-molybdenum  industrial 
complex.  The  charges  were  to  be  placed  at  two  elevations,  2705  m  and  2609  m.  A  general 
scheme  of  the  benches  where  explosions  were  to  be  detonated  and  the  locations  of  charges 
and  the  high-speed  camera  are  shown  at  Hgure  62.  Hgure  63  di^lays  a  rtuq>  of  die  deto¬ 
nation  site,  cross-sections  of  the  bench  and  the  camera  locaticxi. 

Explosive  charges  located  in  73  drilled  holes  were  to  be  aligned  in  S-6  rows  and 
these  charges  were  proposed  to  detonate  within  the  upper  bench  at  the  2705  m  elevations 
(Figure  64).  The  cli^e  in  each  bordiole  were  to  weigh  tqpproximately  350  kg,  each  \xxe^ 
hole  depth  was  14  m,  and  the  total  wei^t  of  the  charges  was  to  be  30  tons.  Ejqilosive 
charges  to  be  located  in  the  38  drilled  holes  were  aligned  in  4  rows  and  these  charges  were 
proposed  to  detonate  within  the  lower  bench  at  the  2690  m  elevation.  Eadi  borehole 
depth  and  each  bcnehole  charge  were  to  be  similar  to  those  on  the  upper  bench.  Total 
weight  of  the  explosive  charges  on  the  lower  bench  was  to  be  15  tons  (Hgure  65). 
According  to  this  project  plan,  the  detonation  time  delay  between  ndghboring  rows  was  to 
be  20  msec. 

The  blasted  block  of  rock  consists  of  marble.  There  are  no  ore  rocks  within  the 
block.  This  rock  category  is  VI  by  explosion  crushing  standards  and  VII  for  drilling.  The 
predicted  distance  for  a  maximum  free  flight  was  3(X)-4(X)  m  and  the  zone  of  seismic  haz¬ 
ard  was  60  m.  The  blast  was  planned  to  be  detonated  on  August  22,  1933  at  6:(X)  a.m 
approximately,  universal  time  (9:(X)  a.m.,  Moscow  time). 
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Hgtire62:  Scheme  of  the  locations  of  the  chaiges  and  camera.  1- suggested  point 
for  filming,  2  -  oversized  pieces  of  rocks  on  the  upper  bench,  3  - 
charge  for  the  upper  bench  at  2705  m,  4  -  oversized  pieces  rocks  on  the 
lower  bench,  5  -  charge  for  the  lower  bench  at  2690  m,  6  -  oversized 
pieces  of  rocks  on  the  lower  bench,  7  -  oversized  pieces  of  rock  on  the 
upper  bench,  8  -  selected  point  for  filming. 
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Figure  6S:  Scheme  of  the  charge  detonations  on  the  lower  bench,  according  to  infor¬ 
mation  obtained  frorn  men  responsible  for  the  explosion  detonations:  1  - 
DSHE-12  blasting  cord,  2  -  electric  detonator. 
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Blast  conditions 


Weadier  conditions  during  the  experiment  on  August  22  were  not  favtvable  widi 
fog  and  low  clouds;  the  air  temperature  was  approximately  4°C  It  had  been  snowing  dur> 
ing  the  night  before  the  experiment  It  took  approximately  2  hours  to  prq)are  for  die  seis> 
mic  measurements.  A  set  of  underground  explosicMis  were  detonated  on  the  same  day 
beftwe  our  experiment  at  9h  9m  39  s  (Moscow  dme).  According  to  the  projea  plan,  die 
total  weight  of  detonated  explosive  charges  was  appmximately  20  tons. 

A  position  for  the  camera  was  selected  only  on  the  explosion  day.  The  new  posi> 
tion  for  the  can^ra  was  chosen  because  the  previously  proposed  position  had  to  be 
rejected  due  to  potential  hazardous  debris.  The  final  camera  position  was  chosen  at  10  h 
07  m  Moscow  time.  The  explosion  was  detonated  at  10  h  27  m  25  sec  (Moscow  time). 
Thus,  we  had  little  time  to  install  the  camera  and  prepare  for  filming.  The  distance  firom 
the  camera  to  the  blast  was  approximately  350  m.  Due  to  lack  of  time,  the  filming  planned 
before  the  experiment  was  not  conduct^  The  camera  was  placed  on  a  tripod  and  was 
switched  on  by  remote  control  through  the  timing  and  device  control  block.  The  camera 
was  turned  on  2  s  before  the  explosion  time. 

After  the  experiment,  we  were  told  by  men  responsible  fix*  the  explosion  detona¬ 
tions  that  the  project  was  changed  during  preparatirms  of  the  charges  within  boreholes. 
The  changes  were  mainly  made  on  the  upper  bench  where  only  5  rows  of  boreholes  were 
detonated  (Hgure  66). 

Meanwhile,  oversized  pieces  of  rocks  were  also  blown  up  on  the  upper  and  lower 
benches.  Figure  62  shows  the  locations  of  the  oversized  rock  pieces.  The  total  weight  of 
the  detonated  explosives  was  only  25  tons  instead  of  the  45  tons  mentioned  in  the  project 
plan,  and  the  time  delays  between  the  drilled  hole  rows  was  40  msec  instead  of  25  msec 
allowed  by  the  project  plan.  According  to  the  specialists  respmisible  for  the  explosion  det¬ 
onations,  the  increased  time  delay  was  caused  by  the  large  dqrtfa  of  the  drilled  holes.  The 
blasting  drcuit  was  made  by  DSHE-12  blasting  cOTd.  The  blasting  cord  was  detonated  by 
an  electric  detonator  using  a  special  blasting  device. 

As  a  result,  we  can  assume  that  changes  in  tire  charges  and  detonation  times  were 
carried  out  arbitrarily,  with  the  deviations  from  the  project  plan  resulting  from  local  pro¬ 
duction  intoests. 

Measurement  Results 

Analysis  of  the  developed  film  showed  that  the  camera  had  a  short  break,  approxi¬ 
mately  120  msec,  immediately  after  origin  time.  Tbis  break  was  probably  caused  by  seis¬ 
mic  vibrations  generated  by  the  explosion  on  the  lower  bench.  For  this  reason,  explosions 
detonated  in  the  upper  and  lower  benches  were  not  filmed  during  the  first  120  msec  after 
the  origin  time.  The  distance  from  the  camera  to  the  explosions  on  the  upper  bench  was 
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Figure  66:  Scheme  ol  ilic  charge  dcioiiaiions  on  the  upper  bench,  according  to  infor¬ 
mation  obtained  from  men  responsible  for  the  explosion  detonations 
(locations  of  detonated  drill  holes  shown  by  the  solid  lines):  1  -  DSHE-12 
blasting  cord,  2  -  electric  detonator. 
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approximately  350  m  (geodedc  measurements  were  not  made).  The  camera  focus  was  35 
mm  and  the  filming  speed  was  70  exposure  per  second.  The  time  between  exposures  was 
14.2  sec  and  signals  with  a  frequency  of  50  Hz  were  used  for  time  calibration.  At  the  end 
of  this  report  are  photos  of  the  explosion  evolution.  Time  measurements  shown  on  these 
photos  were  elapsed  times  after  from  the  origin  time  of  the  explosion  on  the  uf^rer  bei^ 

Analysis  of  the  film  shows  detonations  with  chaotic  and  inaccurate  time  delays. 
This  phenomenon  is  probably  caused  by  variations  in  length  of  the  blasting  cmd.  IMs  will 
consider  the  explodon  evolution  in  det^. 

One  can  see  the  evolution  of  explosions  detonated  within  the  first  row  of  the  upper 
bench  on  the  photo  at  4-30  msec  A  white  cloud,  generated  by  blowing  up  tire  overazed 
pieces  of  rocks,  obscures  some  of  the  nearest  drill  holes. 

The  photo  made  approximately  42  msec  after  the  origin  time  shows  explosions 
detonated  within  the  second  row  of  the  drill  holes.  The  following  photos  show  the  explo¬ 
sion  evolution  vs.  time.  The  photos  made  at  approximately  56  msec  and  70  msec  show  the 
evolution  of  the  explosion  detonated  within  the  second  row  of  drill  holes.  Photos  at  84 
msec  and  98  msec  show  the  detonation  of  the  explosions  within  the  third  row  of  the  drill 
holes.  At  1 12  msec  the  photo  shows  the  evolution  of  the  explosions  detonated  within  tiie 
third  and  fourth  rows  of  the  drill  holes  and  the  beginning  of  the  detonation  of  the  oversized 
pieces  of  rocks.  The  photo  at  126  msec  shows  the  detonation  within  the  fourth  row,  the 
explosion  evolution  for  the  oversized  pieces  of  rock,  and  the  beginning  of  the  explosion 
detonation  within  the  fifth  row  of  the  drill  holes.  At  140  msec  the  photo  shows  the  whole 
blast  evolution.  The  initiM  velocity  of  tiie  cloud  evolution  is  400  n^s. 

It  is  impossible  to  retrieve  the  actual  scheme  of  blasting  because  the  film  analysis 
shows  that  the  designed  plan  was  not  followed  during  the  experiment  We  can  assume  that 
the  blasting  plan  was  changed  during  the  test  in  the  tow  locations.  Moreover,  blasting  of 
the  oversized  pieces  of  rock  created  additional  difficulties  for  film  analysis  of  the  sequence 
of  explotions  in  drill  holes.  The  camera  location  should  be  chosen  based  on  the  exact 
blasting  scheme  and  explosion  features,  in  particular,  for  the  overazed  pieces  of  rocks. 

Analysis  of  the  film  results  showed  that  the  camera  should  be  attached  to  a  firmer 
foundation  such  as  a  buried  wooden  bar.  Two  or  three  cameras  filming  fiom  different 
directions  and  having  a  greater  speed  should  be  used  to  measure  more  accurately  tiie  time 
delays  of  the  explosions. 

The  explosion  improved  visibility,  probably  because  of  fog  condensation  caused 
by  the  influence  of  the  shock  wave.  An  explosion  detonated  in  a  drill  hole  in  the  forth  row 
was  observed  on  the  photo  made  at  84  msec  after  the  origin  time;  this  confirms  our  conclu¬ 
sion  on  the  chaotic  character  of  the  blasting. 
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Conclusions 


1.  The  explosion  filming  confinns  the  chaotic  sequence  of  tiie  detonations  in  the 
drill  holes  and  deviations  of  the  actual  blasting  scheme  fixxn  the  design  plan. 

2.  The  camera  location  should  be  chosen  firom  the  exact  blasting  sdieme  and  die 
explosion  features.  Geodetic  control  of  tiie  camera  location  and  exact  measurements  of 
the  distance  firom  die  camera  to  the  explosion  are  necessary  for  more  accurate  analysis  of 
die  explorion  filming  results. 

3.  Several  cameras  filming  firom  different  directions  and  having  a  greater  qieed 
should  be  used  to  measure  the  explosion  time  delays  more  accurately. 
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